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Time-dependent density functional theory (TDDFT) is presently enjoying enormous popularity 
in quantum chemistry, as a useful tool for extracting electronic excited state energies. This article 
explains what TDDFT is, and how it differs from ground-state DFT. We show the basic formalism, 
and illustrate with simple examples. We discuss its implementation and possible sources of error. 
We discuss many of the major successes and challenges of the theory, including weak fields, strong 
fields, continuum states, double excitations, charge transfer, high harmonic generation, multipho- 
ton ionization, electronic quantum control, van der Waals interactions, transport through single 
molecules, currents, quantum defects, and, elastic electron-atom scattering. 
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I. INTRODUCTION 

Ground-state density functional theory [H,[l,[3[ has be- 
come the method of choice for calculating ground-state 
properties of large molecules, because it replaces the in- 
teracting many-electron problem with an effective single- 
particle problem that can be solved much more quickly. 
It is based on rigorous theorems [H H, E| and a hierar- 
chy of increasingly accurate approximations, such as the 
local density approximation (LDA) , generalized gradient 
approximations (GGA's) [1110], an d hybrids of exact ex- 
change with GGA0]. For example, a recent ground-state 
calculation [i| for crambin (C21Q2H317N55O64S6), a small 
protein, using TURBOMOLEpLCj on a 1.5 GHZ HP ita- 
nium workstation took just 6h52m, extraordinarily fast 
for 2528 electrons. But, formally, ground-state density 
functional theory predicts only ground-state properties, 
not electronic excitations. 

On the other hand, time-dependent density functional 
theory (TDDFT) 0, El El EI El applies the same 
philosophy to time-dependent problems. We replace 
the complicated many-body time-dependent Schrodingcr 
equation by a set of time-dependent single-particle equa- 
tions whose orbitals yield the same time-dependent 
density. We can do this because the Runge-Gross 
theorem [l6j proves that, for a given initial wavefunc- 
tion, particle statistics and interaction, a given time- 
dependent density can arise from at most one time- 
dependent external potential. This means that the time- 
dependent potential (and all other properties) is a func- 
tional of the time-dependent density. 

Armed with a formal theorem, we can then define time- 
dependent Kohn-Sham (TDKS) equations that describe 
non-interacting electrons that evolve in a time-dependent 
Kohn-Sham potential, but produce the same density as 
that of the interacting system of interest. Thus, just as in 
the ground-state case, the demanding interacting time- 
dependent Schrodinger equation is replaced by a much 
simpler set of equations to propagate. The price of this 
enormous simplification is that the exchange-correlation 
piece of the Kohn-Sham potential has to be approxi- 
mated. 

The most common time-dependent perturbation is a 
long-wavelength electric field, oscillating with frequency 
lu. In the usual situation, this field is a weak perturbation 
on the molecule, and one can perform a linear response 
analysis. From this, we can extract the optical absorption 
spectrum of the molecule due to electronic excitations. 
Thus linear response TDDFT predicts the transition fre- 
quencies to electronic excited states and many other 
properties. This has been the primary use of TDDFT 
so far, with lots of applications to large molecules. 

Figure [1] compares TDDFT and experiment for the 
electronic CD spectrum of the chiral fullerene C76. A 
total of 240 optically allowed transitions were required 
to simulate the spectrum. The accuracy is clearly 
good enough to assign the absolute configuration of Cjq. 
TDDFT calculations of this size typically take less than 
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FIG. 1: TDDFT calculation and experiment for the electronic 
CD spectrum of fullerene (y4)-C76. TDDFT calculations were 
performed with the BP86 functional and an augmented SVP 
basis set [l7l |. The RI-J approximation together with TZVP 
auxiliary basis set s [ lal was used. Experimental data (in 
CH2CI2) are from QJp 



a day on low-end personal computers. 

A random walk through some recent papers using 
TDDFT gives some feeling for the breadth of applica- 
tions. Most are in the linear response regime. In inor- 
ganic chemistry, the optical response of man y tr ansition 
metal cc^lexes^jM [M |1 H, H H M H H, 
I30L l3ll [32l [33l . 134 |35| has been calculated, and even 
some X-ray absor p tion 13a . 13711 . In organic chemistry, 
heterocycleslH [M, H| Si El El among others [3, El 
EH have been examined. Other ex amp les include the 
response of thiouracilpTl], s-tetrazine|48l|. and annulated 
porphyrins [49j. We also see TDDFT's use in studying 
various fullerenes [M HH, [M l53l. 154 
TDDFT is finding many uses [561151 

0, [M HU . DNA bases are under examination, and an 
overview of their study may be found in Ref . (6?J • In pho- 



. In biochemistry 
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5J,|5J,|60|,|61|,|6 



tobiology, potential energy curves for the trans-cis photo- 
isomerization of protonatcd Schiff base of retinal[68j have 
been calculated. Large calculations for green and blue 
fluorescent proteins have also been performed [M Ffol ] . 
Doing photochemistry with TDDFT [71(, properties of 
chromophoresEl El El, El, El and dyesfzl El El, 
8(1 l8ll . [H, |83| have been computed. For these and 
other systems, there isgreat interest in charge-transfer 
excitations 0, [M HE H3, H, H IM HH, HI , but (as we 
later discuss) intermolecular charge transfer is a demand- 
ing problem for TDDFT. 

Another major area of application is clusters, large 
and small, covalent and metallic, and every t hing 

96IJ97I, m m DM Em [102 
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inbetween 

ITo5 do 



[9I Hi 



eludin g M e t-Ca rs 111 

sBS [Til, 



effects 



. Se v eral s tudie s inc l ude solven t 

[113, El EM EM EIH, Ell, 

one exampl e be ing the behavior of metal ions in ex- 
plicit water jl23j . TDDFT in linear re s pons e can also 
be used to examine chiralityfUl EM EM EH; in- 
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eluding calculatin g bo th electric and magnetic circular 
dichroismdl fl28l . 1l29l Il30l. [l3lL [l3^ . and has been ap- 
plied to both helical arom atics 133j and to artemisinin 
complexes in soluti on|l34l. There have also b een a ppli- 
cations in materials jl35 , ll3d | and quantum dots jl37l | but, 
as discussed below, the optical respon se of bulk solids re- 
quires some non-local approximations [l38| . 



TABLE I: Table of acronyms and abbreviations. 



Beyond the linear regime, there is also growing inter est 
in second- and third-order response [l39l EM EH EH 

m 

all these fields. In particular t he fie l d of non-l inear optics 
has been heavily investigated jl43l . 1 144 , Il45ll. especially 
the phenomenon of two phot on absorption [l46l Il4q . Il47l 

EH, Eli, EH EIH, E13, EH- 

In fact, TDDFT yields predictions for a huge variety 
of phenomena, that can largely be classified into three 
groups: (i) the non-pcrturbative regime, with systems 
in laser fields so intense that perturbation theory fails, 

(ii) the linear (and higher-order) regime, which yields 
the usual optical response and electronic transitions, and 

(iii) back to the ground-state, where the fluctuation- 
dissipation theorem produces ground-state approxima- 
tions from TDDFT treatments of excitations. 



AIT 

ALDA 
A 
B88 
B3LYP 
CASPT2 
CC 
CIS 
cc 
cxt 
EXX 
GGA 
HK 
H 

HXC 
KS 
LB94 
LDA 
LSDA 
LHF 
LYP 
MAE 
OEP 
PBE 
PBEO 
RG 
RPA 
SPA 
TDKS 
XC 



Asymptotically corrected 
Adiabatic LDA 
Adiabatic 
Becke GGA 

Hybrid functional using Bcckc exchange and LYP correlation 
Complete active space 2 nd order perturbation theory 
Coupled cluster 

Configuration-interaction singlets 

electron-electron 

external 

Exact exchange 

Generalized gradient approximation 

Hohcnbcrg-Kohn 

Hartree 

Hartrcc plus exchange-correlation 
Kohn-Sham 

van Leeuwen-Baerends asymptotically corrected functional 
Local density approximation 
Local spin density approximation 

Localized Hartrcc-Fock (accurate approximation to EXX) 

Lee- Yang-Parr correlation 

Mean absolute error 

Optimised effective potential 

Pcrdcw-Burkc-Ernzcrhof GGA 

Hybrid based on PBE 

Rungc-Gross 

Random phase approximation 
Single pole approximation 
Time-dcpcndcnt Kohn-Sham 
Exchange-correlation 



A. Overview 

This work focuses primarily on the linear response 
regime. Throughout, wc emphasize the difference be- 
tween small molecules (atoms, diatomics, etc.) and the 
larger molecules that are of greater practical interest, 
where TDDFT is often the only practical first-principles 
method. We use napthalene (CioHg) as an example 
to show how the selection of the basis set and of the 
XC functional affects excitation energies and oscillator 
strengths computed using TDDFT. Small molecules are 
somewhat exceptional because they usually exhibit high 
symmetry which prevents strong mixing of the KS states 
due to configuration interaction; also, basis set require- 
ments are often exacerbated for small systems. Naphtha- 
lene is large enough to avoid these effects, yet reasonably 
accurate gas phase experiments and correlated wavefunc- 
tion calculations are still available. 

We use atomic units throughout (e 2 = h = m e = 1), 
so that all energies are in Hartrees (1 H~27.2 eV~627.5 
kcal/mol) and distances in Bohr (~0.529 A) unless oth- 
erwise noted. For brevity, we drop comma's between ar- 
guments wherever the meaning is clear. In DFT and 
TDDFT, there is a confusing wealth of acronyms and 
abbreviations. TableUis designed to aid the readers nav- 
igation through this maze. 

The content of this review is organized as follows. Sec- 
tions [ill and IIIII cover the basic formalism of the the- 
ory, that is needed to understand where it comes from, 
why it works, and where it can be expected to fail. Sec- 
tion [IV] is all about details of implementation, especially 
basis-set selection. On the other hand, section [V] is de- 
voted to performance, and analyzing the sources of er- 
ror in the basis-set limit. In section IVI1 we then look 



at a few atoms in microscopic detail: this is because 
we know the exact ground-state Kohn-Sham potential 
in such cases, and so can analyze TDDFT performance 
in great depth. Section IVTT1 is devoted to the many at- 
tempts to go beyond standard functional approximations, 
and especially discusses where such attempts are needed. 
The last substantial section, section I Villi covers topics 
outside the usual linear response approach to excitations, 
including ground-state functionals derived from TDDFT, 
challenges for strong fields, and transport through single 
molecules. Section [IXI is a summary. 



II. GROUND-STATE REVIEW 

In this section, we review ground state DFT rather 
quickl y. F or a more comprehensive review, we recom- 
mend |l54j . Many of the results discussed here are re- 
ferred to in later sections. 



A. Formalism 

Ground-state DFT is a completely different approach 
to solving the many-electron problem than the traditional 
solution of the Schrodinger equation. The Hohenberg- 
Kohn (HK) theorem Q of 1964 states that for a given non- 
degenerate ground-state density n(r) of Fermions with a 
given interaction, the external potential v CK t (r) that pro- 
duced it is unique (up to an additive constant). Hence 
if the density is known, then v cxt (r) is known and so H, 
the Hamiltonian, is known. From this and the number of 
particles (determined by the integral of the density), all 
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properties of the system may be determined. In partic- 
ular, the ground-state energy of the system E would be 
known. This is what we mean when we say these prop- 
erties are functionals of the density, e.g., E[n]. It was 
later shown that this holds even for degenerate ground- 
states and modern DFT calculations use an analo- 
gous theorem applied to the spin densities, n a (r), n()(r). 
Where a, [3 = ±4 respectfully. 

The total energy for TV electrons consists of three parts: 
the kinetic energy T[$f], the electron-electron interaction 
V cc [\&] , and the external potential energy V ext [&} , each of 
which is defined below: 



i=l 
N N 



N 



V«t[*] = (*|X>«t(rO|*) 

i 

By the Rayleigh-Ritz principle: 



E = min(*| J ff|^) 

= min(T[*] + ^e[*] + Kxt[*]) 



(1) 

(2) 
(3) 



(4) 



If we s imp ly rewrite the minimization as a two step 
process [155} : 



£= min min (T[tf] + VJp] + V ext fi>}) 



where the inner search is over all interacting wavefunc- 
tions yielding spin densities n a ,np. We may pull the last 
term out of the inner minimization: 

E = min \F[n a ,np]+y I d 3 r u 0X t CT (r) n CT (r) 



min (E[n a ,np\) 



(5) 



where 



F[n ai n p ] = min (T[tf] + V ee [V}) (6) 

= Tfna,^] + Vee[n a ,np] , (7) 

is a universal functional independent of f ex to-(r). 

Minimizing the total energy density functional, Eq. 
(JSJ), for both spin densities by taking the functional 
derivative 6/Sn a , and using the Euler-Lagrangc multi- 
plier technique leads to the equation: 



6F[n a ,np] 
5n a 



Wcxta-(r) = fl , 



(8) 



where /i is the chemical potential of the system. 

Next we imagine a system of non-interacting electrons 
with the same spin densities. Applying the HK theo- 
rem to this non-interacting system, the potentials, i> S(T (r), 
that give densities n CT (r) as the ground-state spin den- 
sities for this system are unique. This is the fictitious 
Kohn-Sham (KS) systemQ, and the fully interacting 
problem is mapped to a non-interacting one which gives 
the exact same density. Solving the KS equations, which 
is computationally simple (at least compared to the fully 
interacting problem which becomes intractable for large 
particle numbers), then yields the ground-state density. 
The KS equations are 



"2 V +^( r ) ) <M r ) 



</W r ) 



with spin densites 



(r)=£|^ CT (r)| 5 

3=1 



(9) 



(10) 



where v sa ,v s /3 are the KS potentials and N a is the 
number of spin a electrons , (N a + Np = N). 

In Fig [2J we plot the exact density for the He atom 
from a highly accurate wavefunction calculation, and 
below we plot the exact KS potential for this system. 
One can see that the KS potential is very different from 
the external potential. This is due to the fact that 
the KS single effective potential for the non-interacting 
system must give the correct interacting electron density. 
Because the coulomb repulsion between the electrons 
shields the nucleus, and makes the charge density decay 
less rapidly than e~ 4r , the KS potential is shallower 
than i>ext(r)- 

To derive an expression for v sa (r) we note that the 
Euler equation that yields the KS equations is: 



8T s [n a ,np\ 



+ v sa (r) = fi 



(11) 



Here T s is the kinetic energy of the KS electrons, 

r s = EE / <*V i|VM*)l a ■ 

If we rewrite F[n a ,np\ in terms of the KS system: 

F[n a ,nf3]=T s [n a ,np] + U[n]+E xc [n a ,n[3] , (12) 
where U[n] is the Hartrec energy, given by 

■fi.AJ n CT (r)rv(r') 



r - r' 



(13) 
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and the exchange-correlation (XC) energy is defined by 
Eq. © and Eq. (JT2J: 

E xc [n a ,np] = T[n a ,np]-T s [n a ,np]+V ee [n a ,np]-U[n] ■ 

(14) 

Inserting this into Eq. ([8]) and comparing to Eq. (fTTj) 
gives a definition of the KS potential: 



v s<J (r) = Uext(r) + u H (r) + v XCc7 (r) 



(15) 



where the Hartree potential is the functional derivative 
of U[n}: 



u H (r) 



fl/[n] 
Sn(r) 



E 



<*V ?^ , (16) 



while the XC potential is given by the functional deriva- 
tive of the XC energy: 



8E[n a ,np] 
Sn a (r) 



(17) 



This then closes the relationship between the KS sys- 
tem and the original physical problem. Once E xc [n a , np] 
is known exactly or approximated, v xca (r) is determined 
by differentiation. The KS equations can be solved 
self-consistently for the spin densities and orbitals, and 
the total energy found by inserting these into the total 
energy functional E = T s + U + E xc + V eX f Unfortu- 
nately E xc [n a ,np] is not known exactly and must be 
approximated. There exists a functional soup of many 
different approximations of varying accuracy and compu- 
tational cost. Many of these are discussed in Section lll Bl 

In Fig [2j when the KS equation is solved with this 
exact potential, the HOMO level is at -24.592 eV. This 
is minus the ionization energy for Helium. In exact 
DFT, Koopma n's t heorem, which states / = —£homo, 
is exactly true jl57| . In ground-state DFT, this is the 
only energy level of the fictitious KS system that has an 
immediate physical interpretation. 




FIG. 2: Top panel - exact radi al density for the He atom 
found via the QMC method )l56l ]. Bottom panel - The ex- 
ternal and KS potentials for the He atom. The KS potential 
is found by inverting the KS equations using the exact KS 
orbitals (easily found for He if exact density is known). 



Before leaving the ground-state review, we menti o n th e 
optimized effective potential (OEP) method |l58l . I X 59j ] . 
Here the XC functional is written as a functional of the 
KS orbitals (which in turn arc functionals of the den- 
sity). The exchange energy is then given by the familiar 
HF definition. 



EXX gives useful features such as derivative discon- 
tinuities and the correct asymptotic decay of the KS 
potential jl6C| . As we will see in Section TV B[ these are 
important for TDDFT linear response. 



E x 



1 N 

EE 



drdr' 



</ar)</>Ur')<Mr)<Mr') 



(18) 

However in contrast to HF, a single effective potential 
i>* x (r) is found via the chain rule. Accurate orbital 
dependent functionals for the correlation energy are 
extremely difficult to find, so often only exchange is 
used. In DFT, this is called exact exchange (EXX), 
since exchange is usually only treated approximately. 



B. Approximate Functionals 

In any ground-state DFT calculation, we must use 
approximations for the functional dependence of the 
XC energy on the spin densities. There now exists 
a hierarchy of such approximations. The simplest 
of these is the local density approximation (LDA), 
where the XC energy at a point r' is calculated as if 
it were a uniform electron gas with the spin densities 



6 



n a = n a (r') in a constant positive background. The 
exchange energy for the uniform gas can be deduced 
analytically, but the correlation contribution is found 
using a combination of many-body theory and highly 
accurate Monte C arlo simulations fo r the electron gas of 
different densities [TH GH, EEH Ell ■ 

LDA works remarkably well, given the vast different 
between homogeneous electron gases and atoms or 
molecules. However total energies are generally underes- 
timated. Typically the XC energy is underestimated by 
about 7%. When the performance of LDA is examined 
carefully, this comes about via a nice (but not completely 
accidental) cancellation of errors between the exchange 
part (underestimated by about 10%) and correlation 
(overestimated by 200% — 300%), which is typically 4 
times smaller than exchange. 

An obvious improvement to LDA would be to include 
information about how rapidly the density is changing 
via its gradient. This leads to the generalized gradient 
approximation (GGA). In the original Kohn-Sham 
paper of 1965, the simplest of these was suggested. 
The gradient expansion approximation (GEA) is found 
by ex amin ing the slowly varying limit of the electron 
gas f2i [l65j . However it was soon found that GEA failed 
to improve on the accuracy of LDA and sometimes made 
things worse. It was not until the late 80's that accurate 
GGA's were constructed. The most popular of these are 
BLYP (B880 for exchange and LYPpj for correlation) 
and PBE0- These generally reduce atomization errors 
of LDA by a factor of 2 — 5. 

PBE is a functional designed to improve upon the 
performance of LDA without losing the features of 
LDA which are correct. As such it reduces to LDA 
for the uniform electron gas. A GGA should also 
satisfy as many exact conditions as possible, such as 
the Licb-Oxford bound or the form of the exchange 
energy in the slowly varying limit of the electron gas. 
In this regard, PBE is a non-empirical functional where 
all parameters are determined by exact conditions. 
Because of its ability to treat bulk metals, it is the 
functional of choice in solid-state calculations and 
increasingly so in quantum chemistry. When choosing 
a GGA, using PBE or not PBE should no longer 
be a question. (Although the crambin calculation of 
the introduction used BP86, for reasons explained in 0].) 

Finally, hybrid functionals mix in some fraction of 
exact exchange with a GGA. This is the Hartrcc-Fock ex- 
change integral, Eq. (fT8"|) , evaluated with the KS orbitals 
(which are functionals of the density). Only a small 
fraction of exact exchange (20% — 25%) is mixed in, in 
order to pre serve the cancellation of errors which GGA's 
make use of [l66l | . The most widely used functional in 
chemistry is the hybrid function B3LYP, whi ch contains 
3 experimentally fitted parameters @, H, Il67| (although 



the parameter in B88 has recently been derived [l68l ]L 
Other hybrid functionals include PBE0, where 2 5% o f 
exact exchange is mixed in with the PBE functional [l6^]. 
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FIG. 3: Exact and LDA KS potentials for the He atom. While 
the exact potential falls off as — 1/r, the LDA decays much 
too quickly. This is common for nearly all present functionals 
and has major consequences for TDDFT. 

A less well-known feature to users of ground state DFT 
is that while their favourite approximations yield very 
good energies (and therefore structures, vibrations, ther- 
mochemistry, etc.) and rather good densities, they have 
poorly behaved potentials, at least far from nuclei. Fig- 
ure [3] illustrates this for the He atom, showing the LDA 
potential compared to the exact KS potential. While the 
potential is generally good in the region r < 2, it decays 
much too fast far from the nucleus. The true KS poten- 
tial falls off as — 1/r whereas LDA decays exponentially. 
Hence the KS eigenvalues and eigenvalues will be poor 
for the higher levels. To understand why poor potentials 
do not imply poor energies (and why t hese potentials are 
not as bad as they look), see Ref. [l70j |. But, as we 
shall see in section [V] this has major consequences for 
TDDFT. 

Over the past decade, the technology for treating 
orbital-dependent functionals has deve loped, and such 
functionals help cure this problem [158| |. Th i s is called 
the optimized effective potential (OEP)[l7j [HI Gil]. 
The first useful orbital functional was the s elf-interaction 
corrected LDA of Perdew and Zunger |l62j . More gener- 
ally, the OEP method can handle any orbital-dependent 
functional including treating exchange exactly. Orbital- 
dependent functionals naturally avoid the self-interaction 
error that is common in density functionals. An (al- 
most) exact implementation o f th e OE P equations is 
localized Hartree-Fock (LHF) [l74 Il75fl . available in 
TURBOMOLEpl]. 



C. Basis Sets 

To actually solve the KS equations, the KS orbitals 
4> p<7 {y) are expanded in a finite set of basis functions 
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4> pa {vt) = ^2 G vva x„(r). (19) 

V 

The most common choice by far for the basis functions in 
quantum chem istry are atom-centered contracted Carte- 
sian Gaussians [l76j |. 

X„(r) ^c^M^M^U^-fcfr-R.) 2 . (20) 

i 

l x (v)i lyiv), and l z {v) are positive integers or zero, and 
l{v) = lx(v) + ly{ v ) + 1>z{ v ) is somewhat loosely called l- 
quantum number of %„. (I = 0, 1, 2, 3, . . . corresponds to 
s,p,d, f, . . . type Cartesian Gaussians.) The exponents 
Cip, and the contraction coefficients a v are optimized in 
atomic calculations. Other common basis functions in 
use are Slater type orbitals, plane waves, or piecewise 
defined functions on a numerical grid. 

The approximation of the orbitals </> p(T (r) by a finite 
linear combination of basis functions (also called LCAO, 
linear combination of atomic orbitals), Eq. (fl9l) . leads 
to a finite number of MOs. Thus, the KS equations 
and all derived equations are approximated by finite- 
dimensional matrix equations. These equations can be 
treated by established numerical linear and non-linear al- 
gebra methods. When the basis set size is systematically 
increased, the computed properties converge to their ba- 
sis set limit. 

In a finite basis set, all operators become finite matri- 
ces; the matrix elements are integrals, e.g., 

H Mm [n] = J d\ x ^)HMxA*)- (21) 

The calculation and processing of such integrals is the 
main effort in virtually all DFT calculations. Gaussian 
basis functions have the distinct advantage that most in- 
tegrals can be evaluated analytically, and that they are 
spatially local. The latter implies that many integrals 
vanish and need not be calculated. Whether a certain 
integral vanishes or not can be decid ed in advance by 
so-called pre-screening techniques [l77j . 

To illustrate the effect of choosing various basis sets, 
we show in Fig. [TH the reaction energy for naphthalene 
combustion in the gas phase: 

CioH 8 + 120 2 -> 10CO-2 + 4H 2 (22) 

The basis sets are listed in order of increasing size, 
and are well-know in quantum chemistry (and are 
described in detail in section [TV| . We see that hydrogen 
polarization functions (basis sets ending in P) are 
important, because C-H bonds are broken and O-H 
bonds are formed. Augmentation (aug-) with diffuse 
functions somewhat improves the smaller basis-set 
results, but is not economical in this case. Using the 
resolution of the identity for the Coulomb operator (RI) 
saves computational time, with no loss of accuracy. 



TABLE II: Single-point calculations using PBE functional 
for the reaction energy for naphthalene combustion using 
PBE/TZVP/RI geometries. Reference value co mpu ted us- 
ing standard enthalpies of formation(from NIST 178]) using 
thermal and ZVPE corrections at the PBE/TZVP/RI level. 



Basis set 


Negative reaction energy (kcal/mol) 


SV 


916.8 


6-31G* 


i nfin n 
1047.1 


SVP 


1108.6 


aug-SV(P) 
TZVP 


1115.5 
1124.5 


TZVPP 


1131.2 


cc-pVTZ 


1129.0 


aug-TZVP 


1130.2 


aug-TZVP/RI 
QZVP 


1130.2 
1140.3 


Reference Value 


1216.3 



Reasonable results are found with SVP, but convergence 
improves all the way to TZVP. We can see after the 
TZVPP result, the basis set error is below the functional 
error and the result is effectively converged. We have 
reached the stage where adding more orbitals, which 
increases the computational cost, is no longer going to 
drastically improve the result. (On the other hand, the 
crambin calculation mentioned in the introduction is 
very large and so only an SV(P) basis set could be used). 

On first impression, comparison to the reference value 
indicates quite a large error, AE= 76 kcal/mol. However, 
given that 48 electron pair bonds are broken and formed, 
the error per carbon atom, 7.6 kcal/mol, is typical for 
this functional. 



III. TIME-DEPENDENT THEORY 

In this section, we introduce all the basis elements of 
TDDFT, and how it differs from the ground-state case. 

A. Runge-Gross theorem 

The analog of the Hohenbcrg-Kohn theorem for time- 
dependent problems is the Runge-Gross theorem [ill, 
which we state here. Consider N non-relativistic elec- 
trons, mutually interacting via the Coulomb repulsion, in 
a time-dependent external potential. The Runge-Gross 
theorem states that the densities n(rt) and n'(rt) evolv- 
ing from a common initial state ^/q — ^/(t — 0) under the 
influence of two external potentials v ex t(rt) and v' Q7Ct (rt) 
(both Taylor expandable about the initial time 0) are al- 
ways different provided that the potentials differ by more 
than a purely time-dependent (r-independent) function: 

Av oxt (rt) ? c(t) , (23) 
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where 



Av cxt (rt) = Vext(rt) - 



; ext 



(24) 



Thus there is a one-to-one mapping between densities and 
potentials, and we say that the time-dependent potential 
is a functional of the time-dependent density (and the 
initial state). 

The theorem was proven in two distinct parts. In the 
first (RGI), one shows that the corresponding current 
densities differ. The current density is given by 



j(rt) - j(r) mm 



where 



j(r) 



1 N 



r.O+^r-r^V,-) 



(25) 



(26) 



is the current density operator. The equation of motion 
for the difference of the two current densities gives jig: 



dAj(rt) 



Of 



= -n (r)VA« ext (r,0) 



(27) 



If the Taylor-expansion about t = of the difference of 
the two potentials is not spatially uniform for some order, 
then the Taylor-expansion of the current density differ- 
ence will be non-zero at a finite order. This establishes 
that the external potential is a functional of the current 
density, v ext [ j, *o](r, *)• 

In the second part of the theorem (RGII), continuity 
is used: 



dn(rt) 

at 



-V-j(rt) 



(28) 



which leads to: 

d 2 An{rt) 



dt 2 



V- (n o (r)VA« ext (r,0)) (29) 



t=o 



Now, suppose At! ext (r, 0) is not uniform everywhere. 
Might not the left-hand-side still vanish? A ppare ntly not, 
for real systems, because it is easy to show jl79l |: 



d 3 r Aw ext (r,0)V • (n o (r)VA Wext (r,0)) 



[V • (Aw oxt (r,0)no(r)VA% x t(r,0)) 



-TiolVAuextCr.O)! 5 



(30) 



Using Green's theorem, the first term on the right van- 
ishes for physically realistic potentials (i. e., potentials 
arising from normalizable external charge densities), 
because for such potentials, Au ex t(r) falls off at least as 
1/r. But the second term is definitely negative, so if 
Aii ex t(r,0) is non- uniform, the integral must be finite, 
causing the densities to differ in 2nd order in t. This 
argument applies to each order and the densities n(r,t) 



and n'(r,t) will become different infinitcsimally later 
than t. Thus, by imposing this boundary conditions, we 
have shown that v ex t[n, ^o](jt). 

Notes: 

• The difference between n(rt) and n'(rt) is non- 
vanishing already in first order of Av ex ± (rt) , ensur- 
ing the invcrtibility of the linear response operators 
of section Ull CI 

• Since the density determines the potential up to 
a time-dependent constant, the wavefunction is in 
turn determined up to a time-dependent phase, 
which cancels out of the expectation value of any 
operator. 



• We write 



■Uext[n;*o](rt) 



because it depends on both the history of the den- 
sity and the initial wavefunction. This functional is 
a very complex one, much more so than the ground- 
state case. Knowledge of it implies solution of all 
time-dependent Coulomb interacting problems. 

• If we always begin in a non-degenerate ground- 
state [HO, EHI, the initial-state dependence can be 
subsumed by the Hohenberg-Kohn theorem[l[, and 
then fext(rf) is a functional of n(rf) alone: 

• A spin-dependent generalization exists, so that 
Vext(rt ) wi ll be a functional of the spin densities 
n a , n^[l82|. 



B. Kohn-Sham equations 

Once we have a proof that the potential is a functional 
of the time-dependent density, it is simple to write the 
TD Kohn-Sham (TDKS) equations as 

i^^=(~+tU»](rt))Mrt) , (31) 

whose potential is uniquely chosen (via the RG theorem) 
to reproduce the exact spin densities: 



r (rtH£|<Mrt)| 2 , 

3=1 



(32) 



of the interacting system. We define the exchange- 
correlation potential via 

/n(r't) 
d 3 r' -A_ L +Vxca (rt) , (33) 

where the second term is the familiar Hartree potential. 

Notes: 
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• The exchange-correlation potential, v XC(J (rt), is in 
general a functional of the entire history of the 
densities, n a {rt), the initial interacting wavefunc- 
tion ^(0) , an d the initial Kohn-Sham wavefunc- 
tion, 'l>il)i;l >l . But if both the KS and interacting 
initial wavefunctions are non-degenerate ground- 
states, it becomes a simple functional of n a (rt) 
alone. 

• By inverting the single doubly-occupied KS equa- 
tion for a spin-unpolarized two-electron system, it 
is quite straightforward (but technically demand- 
ing) to find the TDKS potential from an exact 
time- dependent den sity, and has been done several 
timesfl8l[T83.[l8l. 



• In practical calculations, some approximation is 
used for v xc (rt) as a functional of the density, and 
so modifications of traditional TDSE schemes are 
needed for the propagation [18f|. 

• Unlike the ground-state case, there is no self- 
consistency, merely forward propagation in time, 
of a density dependent Hamiltonian. 

• Again, contrary to the ground-state, there is no 
central role played by a single-number functional, 
such as the ground-state energy. In fact, an action 
was written down in the RG paper, but extrem- 
izing it was l ater shown not to yield the TDKS 
equations jl87j . 



C. Linear response 

The most common application is the response to a 
weak long-wavevlcngth optical field, 



Sv ext {rt) = -£exp(iwi)z 



(34) 



In the general case of the response of the ground-state to 
an arbitrary weak external field, the system's first-order 
response is characterized by the non-local susceptibility 

S n<7 (rt)=J2 J dt' J d 3 r'x**>{n }(r,r';t-t r )6v cx t*>(r't'). 

(35) 

This susceptibility ^ is a functional of the ground-state 
density, rio(r). A similar equation describes the density 
response in the KS system: 

5n a {rt)=Yl J dt ' J d 3 r' Xs *o'[no}(r,r';t-t') fcw(rV). 

(36) 

Here Xs is the Kohn-Sham response function, constructed 
from KS energies and orbitals: 



<(rr w) = <W > < j— : q — r— 



where q is a double index, representing a transition 
from occupied KS orbital i to unoccupied KS orbital 
a, uj q<J = e a * - e icr , and $ g<r (r) = ^(^^(r). 0+ 
means the limit as 0+ goes to zero from above (i.e., 
along the positive real axis). Thus Xs is completely 
determined by the ground-state KS potential. It is the 
susceptibility of the non-interacting electrons sitting in 
the KS ground-state potential. 

To relate the KS response to the true response, we 
examine how the KS potential in Eq. (|3"3"]) changes: 



5v sa (r,t) = Sv cxta (r, t) + 8v uxca {rt) 



(38) 



Since 8v nxC a (ri) is due to an infinitesimal change in 
the density, it may be written in terms of its functional 
derivative, i.e., 

Sv HXCU (rt) = J2 J rfV / dt' f nxca(T ,(rr',t-t') 5n a ,{v't r ) 

(39) 



where 



.[n ](rr',i-i') 



Sv H 



(ri) 



5n a >{T't') 



(40) 



The Hartree contribution is found by differentiating Eq. 
(El: 



f n [n ](rr',t-t') = 



5{t-t') 



(41) 



while the remainder / XOffff i [no](rr', t — t') is known as 
the XC kernel. 

By the definition of the KS potential, 8n a {rt) is the 
same in both Eq. (|35j) and Eq. (j36|) . We can then 
insert Eq. §§§ into Eq. $o§, equate with Eq. (J35) and 
solve for a general relation for any 5n a (rt). After Fourier 
transforming i n tim e, the central equation of TDDFT 
linear response 188] is a Dyson-like equation for the true 
X of the system: 

Xaa'(rr'w) = xsaa' (tt'u) + ^ J d 3 ?'i J d 3 r 2 Xs**! (niw) 



/xc<T 1( r 2 (rir2w) ) x* 2 *'( r 2r'w 



(42) 



(37) 



Notes: 

• The XC kernel is a much simpler quantity than 
v xc *[n\{Yt) , since the kernel is a functional of only 
the ground-state density. 

• The kernel is non-local in both space and time. The 
non-locality in time manifests itself as a frequency 
dependence in the Fourier transform, / XC[Jff )(rr'w). 
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• If /xc is set to zero in Eq. (|42|) . physicists call 
it the Random Phase Approximation (RPA). The 
inclusion of / xc is an exactification of RPA, in the 
same way the inclusion of v xc (j) in ground-state 
DFT was an exactification of Hartree theory. 

• The Hartree kernel is instantaneous, i.e., local in 
time, i.e., has no memory, i.e., given exactly by an 
adiabatic approximation, i.e., is frequency indepen- 
dent. 

• The frequency-dependent kernel is a very sophisti- 
cated object, since its frequency-dependence makes 
the solution of an RPA-type equation yield the 
exact x (including all vertex corrections at every 
higher order term). It defies physical intuition and 
arguments based on the structure of the TDDFT 
equations are at best misleading. If any argument 
cannot be given in terms of many-body quantum 
mechanics, Eq. (|42|) cannot help. 



or 



• The kernel is, in general, complex, with re al an d 
imaginary parts related via Kramers-Kronig [l89j . 



Next, Casida 190] used ancient RPA technology, to 
produce equations in which the poles of \ are found as 
the solution to an eigenvalue problem. The key is to ex- 
pand in the basis of KS transitions. We write 8n a (jct) 



8n a {vu) = J2 (JVM^tt + (43) 



ilqaq' a' {to) = Uj\ a 8 qq <8 aa > + 2^U q(y UJ q > a > K qaq i a > . (48) 



Oscillator strengths f q may be calculated [190J from the 
normalized eigenvectors using 



f qa = i (|x T S- 1 /2^|2 + ijfs-i/ag^ + ^s-i/a^a) , 

(49) 

where 

S qq ' ~ 8 qq > 8 aa i I w q t . 



> 23 



21) 



Continuum 



Exact KS 



Exact TDDFT 



where q = (a, i) if q = (i, a). This representation is used 
to solve Eq. (|36|) self-consistently using Eq. (j39|) . and 
yields two coupled matrix equations [l9l| : 



A B 



B* A* I w V 1 



-1 



8w 
Sv* 



(44) 

whciG Aqaq'fj' — ^qq' Sa-a-' LUq(j-\- Kqfjqf a f ^ Bq a q> a i — K qaq' 'a' i 
Xqo- — Pqai Yq a — Pqa &nd 

K qaq , a ,{oj) = J dr J dr' $ gCT (r) / HX c ff(I ' (rr'w) $*, CT ,(r') , 

(45) 



with 



8v qa (u) = dr $ gcr (r) 8v c ^ t (ru) 



(46) 



At an excitation energy, Sv = and choosing real KS 
orbitals and since (A — B) is positive definite, we get: 



&qcrq'a'{w) 3 q > a i 



uj 2 d„ 



(47) 



where 



n = (A - B) 1 / 2 (A + B)(A - B) 1 / 2 , 



FIG. 4: Transitions for the Helium atom using in ground-state 
DFT on the left, and TDDFT on the right. In both cases, the 
exact functionals have been used. The results for employing 
the exact XC kernel in TDD FT l inear response are known 
from calculations using Ref. |l92|. In each pair of lines on 
the right, the triplet is the lower. 

Figure [4] shows the results of exact DFT calculations 
for the He atom. On the left, we consider just transitions 
between the exact ground-state KS occupied (Is) to un- 
occupied orbitals. These are not the true excitations of 
the system, nor are they supposed to be. However, ap- 
plying TDDFT linear response theory, using the exact 
kernel on the exact orbitals, yields the exact excitations 
of the He atom. Spin-decomposing produces both singlet 
and triplet excitations. 



D. Approximations 

As in the ground-state case, while all the equations 
above are formally exact, a practical TDDFT calculation 
requires an approximation for the unknown XC potential. 
The most common approximation in TDDFT is the adi- 
abatic approximation, in which 



: v TcA n o]( T )\noAr)=n^rt), (50) 
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i.e., the XC potential at any time is simply the ground- 
state XC potential at that instant. This obviously be- 
comes exact for sufficiently slow perturbations in time, 
in which the system always stays in its instantaneous 
ground-state. Most applications, however, are not in 
this slowly varying regime. Nevertheless, results obtained 
within the adiabatic approximation are remarkably accu- 
rate in many cases. 

Any ground-state approximation (LDA, GGA, hybrid) 
automatically provides an adiabatic approximation for 
use in TDDFT. The most famous is the adiabatic lo- 
cal density approximation (ALDA). It employs the func- 
tional form of the static LDA with a time-dependent den- 
sity: 



ALDA r 

J xca 



~[n](ri) =vZ lt (n a {rt),np{rt)) 



de 



unif 



dn a 



dn a dn„ 



An) 
(51) 

Here e^'l 11 (n a , np) is the accurately known exchange- 
correlation energy density of the uniform electron gas of 
spin densities n j , . For the time-dependent exchange- 
correlation kernel of Eq. l|30]) , Eq. (f5T) leads to 

f^[n }(vt,rt') = <$( 3 >(r-r')<J(*-i') „ 

n a =n aa (r) 

(52) 

The time Fourier-transform of the kernel has no 
frequency-dependence at all in any adiabatic approxima- 
tion. Via a Kramers-Kronig relation, this implies that it 
is purely real (l89j | . 

Thus, any TDDFT linear response calculation can be 
considered as occuring in two steps: 

1. An approximate ground-state DFT calculation is 
done, and a self-consistent KS potential found. 
Transitions from occupied to unoccupied KS or- 
bitals provide zero-order approximations to the op- 
tical excitations. 

2. An approximate TDDFT linear response calcula- 
tion is done on the orbitals of the ground-state cal- 
culation. This corrects the KS transitions into the 
true optical transitions. 

In practice both these steps have errors built into them. 



IV. IMPLEMENTATION AND BASIS SETS 

In this section we discuss how TDDFT is implemented 
numerically. TDDFT has the ability to calculate many 
different quantities and different techniques are some- 
times favored for each type. For some purposes, e.g., 
if strong fields arc present, it can be better to propa- 
gate forward in time the KS orbit als u sing a real space 
grid [HI Il94j or with plane waves |l95j . For finite-order 
response, Fourier transforming to frequenc y sp ace with 
localized basis functions may be preferable [19g. Below, 
we discuss in detail how this approach works, emphasiz- 
ing the importance of basis-set convergence. 



A. Density matrix approach 

Instead of using orbitals, we can write the dynamics of 
the TDKS systems in terms of the one-particle density 
matrix 7 <T (rr / i) of the TDKS determinant. j a (rr't) has 
the spectral representation 



N 



(53) 



i.e., the N a TDKS orbitals are the eigenfunctions of ^y a . 
The eigenvalue of all TDKS orbitals, which is their occu- 
pation number, is always 1, which reflects the fact that 
the TDKS system is non-interacting. Equivalently, ^ a 
satisfies the idempotency constraint 



It 



■ (rr'i) = J dxi7 CT (rrit)7 (T (rir'i). 



(54) 



The normalization of the TDKS orbitals implies that the 
trace of j a be N a . 

Using the TDKS equations (|3"Tj) . one finds that the 
time-evolution of 7 CT is governed by the von-Neumann 
equation 



9 , . 
i— 7 At) 



[H a [n](t),lAt)]- 



(55) 



dt 

where H a [n]{rt) = -V 2 /2 + v sa [n]{rt) is the TDKS one- 
particle Hamiltonian. Although 7^ has no direct physical 
meaning, it provides the interacting density and current 
density: The density is simply 

n a {rt) = 7a(rrt), (56) 

and the current density can be obtained from 

j a (rt) = l(v r -V r /h CT (rr't) . (57) 

11 r'=r 

Thus, one can either propagate the TDKS orbitals using 
the TDKS equations (|3ip , or equivalently one can propa- 
gate the TDKS one-particle density matrix 7^ using the 
von-Neumann equation (|55p . subject to the idempotency 
contstraint (|54|) and normalized to N a . 

In practice, it is often preferable to use j a instead of 
the TDKS orbitals. "f a is unique (up to a gauge transfor- 
mation), while the orbitals can be mixed arbitrarily by 
unitary transformations. Both n a and j a are linear in j a , 
while they are quadratic in the orbitals; also, the TDKS 
equations are inhomogeneous in the orbitals due to the 
density dependence of H a , while they are homogeneous 
in 7^. A response theory based on the TDKS density 
matrix is therefore considerably simpler than one based 
on the orbitals. Finally, the use o f 7 CT is computationally 
more efficient than using orbitals [196J. 



B. Basis Sets 

In response theory, the basis functions X^( r ) are usu- 
ally chosen to be time- independent; for strong fields or 
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coupled electron-nuclear dynamics, time-dependent basis 
functions can be more appropriate. 

C. Naphthalene converged 

Table ITTT1 shows the basis set convergence of the first six 
singlet excitation energies of naphthalene computed us- 
ing the PBE XC functional; the corresponding oscillator 
strengths for some of the transitions are also given. Simi- 
lar basis-set convergence studies on small model systems 
should precede applications to large systems. In prac- 
tice, the systems and states of interest, the target accu- 
racy, the methods used, and the computational resources 
available will determine which basis set is appropriate. 

With a model small molecule, we can find the basis-set 
convergence limit of a method. Both excitation energies 
and oscillator strengths are essentially converged within 
the aug-QZVP basis set. QZVP stands for a quad rupl e 
zeta valence basis set with polarization functions [19!|, 
and the prefix aug- denotes additional diffuse functions 
on non- hydrogen atoms, w hich were taken from Dun- 
ning's aug-cc-pVQZ basis set [200| . For C and H, this cor- 
responds to [8s5p4d3/2g] and [4s3p2dl/], respectively, 
where the numbers in brackets denote shells of contracted 
Gaussian type orbitals (CGTOs), as usual. We will take 
the aug-QZVP results as a reference to assess the effect 
of smaller basis sets. 



D. Double zeta basis sets 

The smallest basis in Table IIIII is of split valence (SV) 
or double zeta valence quality [20l|, without polarization 
functions. This basis set consists of two CGTOs per va- 
lence orbital and one per core orbital, i.e. [3s2p] for C 
and [2s] for H. Anoth er popular double zeta valence ba- 
sis set is 6-31G [202j |. The SV basis set can be used to 
obtain a very rough qualitative description of the low- 
est valence excited states only, e.g. 1 1 Bz u and 1 1 £?2ii- 
Higher and diffuse excitations, such as 1 1 A U , are much 
too high in energy or can be missed completely in the 
spectrum. Since unpolarized basis sets also give poor 
results for other properties such as bond lengths and en- 
ergies, their use is generally discuraged nowadays. 

E. Polarization functions 

By adding a single set of polarization functions to non- 
hydrogen atoms, the SV results for valence excitations 
can be considerably improved, at still moderate compu- 
tational cost. The resulting basis set is termed SV(P) 
and consists of [3s2pld] for C and [2s] for 11 201. The 
basis set errors in the first two valence excitation ener- 
gies is reduced by about 50%. There is also a dramatic 
improvement in the oscillator strength of the dipole al- 
lowed transitions. This is expected from the limiting case 



of a single atom, where the first dipole allowed transition 
from a valence shell of I quantum number l v generally 
involves orbitals with Z-quantum number l v + 1. Basis 
sets of SV(P) or similar quality are often the first choice 
for TDDFT applications to large systems, especially if 
only the lowest states are of interest and/or diffuse exci- 
tations are quenched, e.g. due to a pola r environment. 
The popular 6-31G* basis set [202l |203| has essentially 
the same size as SV(P) but performs slightly poorer in 
our example. 

Adding a single set of p type polarization fun ction s 
to hydrogen atoms produces the SVP basis set [20 lj . 
These functions mainly describe C-H a* type excitation 
in molecules which usually occur in the far UV and are 
rarely studied in applications. In our example, going 
from SV(P) to SVP has no significant effect. This may 
be different for molecules containing strongly polarized 
hydrogen-clement or hydrogen bridge bonds. 

Next, aug-SV(P) is an SV(P) basis set augmented by 
a [lslplc?] set of primitive Gaussian s wit h small expo- 
nents (from Dunning's aug-cc-pVDZ [20d| ). often called 
"diffuse functions". As shown in Table IITTI the effect of 
diffuse augmentation is a moderate downshift of less than 
0.1 eV for the first two singlet excitation energies. This 
behavior is typical of lower valence excited states having 
a similar extent as the ground-state. Our example also 
shows that diffuse functions can have a significant effect 
on higher excitations. An extreme case is the 1 1 A U state 
which is an excitation into the lOau orbital having the 
character of a 3s Rydberg state (of the entire molecule). 
The excitation energy of this state is lowered by more 
than 1 eV upon diffuse augmentation. 

While polarization functions arc necessary for a quali- 
tatively correct description of transition dipole moments, 
additional diffuse polarization functions can account for 
radial nodes in the first-order KS orbitals, which fur- 
ther improves computed transition moments and oscilla- 
tor strengths. These benefits have to be contrasted with 
a significant increase of the computational cost: In our 
example, using the aug-SV(P) basis increases the com- 
putation time by about a factor of 4. In molecules with 
more than 30-40 atoms, most excitations of interest are 
valence excitations, and the use of diffuse augmentation 
may become prohibitively expensive because the large ex- 
tent of these functions confounds integral prescreening. 



F. Triple zeta basis sets 

In such cases, triple zeta valence (TZV) basis sets can 
be a better alternative. The TZVP (def-2-TZVP, Ref. 
[2041 ]') basis set corresponds to [5s3p2dlf] on C and [3slp] 
on H. It provides a description of the valence electrons 
that is quite accurate for many purposes when density 
functionals are used. At the same time, there is a sec- 
ond set of polarization functions on non- hydrogen atoms. 
The excitation energies of valence states are essentially 
converged in this basis set, see Table HTT1 However, diffuse 
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TABLE III: Basis set convergence of the first six singlet excitation energies (in eV) and oscillator strengths (length gauge) of 
naphthalene. The basis set acronyms are defined in the text. iVbf denotes the number of Cartesian basis functions, and CPU 
denotes the CPU time (seconds) on a single processor of a 2.4 GHz Opteron Linux workstation. The PBE functional was used 
for both the ground-state and TDDFT cal culat ions. The ground-state structure was optimized at the PBE/TZVP/RI-level. 
Experimental results were taken from Ref. [197| [. 



Basis set 


1 1 i?3ii 


1 1 B 2u (Osc. Str.) 


2 M 9 


1 l B la 


2 1 B Zu (Osc. Str.) 


1 l A u 


N hi 


CPU 


SV 


4.352 


4.246 (0.0517) 


6.084 


5.254 


5.985 (1.1933) 


6.566 


106 


24 


SV(P) 


4.272 


4.132 (0.0461) 


5.974 


5.149 


5.869 (1.1415) 


6.494 


166 


40 




4.290 


A 1 C A 

4.1o4 


O.U21 


o.loo 


0.902 


( .013 


loo 


4U 


SVP 


4.262 


4.125 (0.0466) 


5.960 


5.136 


5.852 (1.1402) 


6.505 


190 


48 


aug-SV(P) 


4.213 


4.056 (0.0417) 


5.793 


4.993 


5.666 (1.1628) 


5.338 


266 


168 


TZVP 


4.209 


4.051 (0.0424) 


5.834 


5.030 


5.715 (1.1455) 


6.215 


408 


408 


TZVPP 


4.208 


4.050 (0.0425) 


5.830 


5.027 


5.711 (1.1464) 


6.231 


480 


568 


cc-pVTZ 


4.222 


4.064 (0.0427) 


5.870 


5.061 


5.747 (1.1355) 


6.062 


470 


528 


aug-TZVP 


4.193 


4.031 (0.0407) 


5.753 


4.957 


5.622 (1.1402) 


5.141 


608 


2000 


aug-TZVP/RI 


4.193 


4.031 (0.0407) 


5.752 


4.957 


5.621 (1.1401) 


5.142 


608 


400 


QZVP 


4.197 


4.036 (0.0416) 


5.788 


4.989 


5.667 (1.1569) 


5.672 


1000 


6104 


aug-QZVP 


4.192 


4.029 (0.0406) 


5.748 


4.954 


5.616 (1.1330) 


5.071 


1350 


28216 


expt. 


3.97, 4.0 


4.45, 4.7 (0.102, 0.109) 


5.50, 5.52 


5.28, 5.22 


5.63, 5.55, 5.89 (1.2, 1.3) 


5.6 







states are too high in energy. There is very little change 
going to the TZVPP basis, which differs from TZVP only 
by an additional set of po lariz ation functions on H. Dun- 
ning's cc-pVTZ basis set j205| performs similar to TZVP 
and TZVPP. However, Dunning basis sets are based on 
a generalized contraction scheme for valence orbitals, as 
opposed to the segmented contracted SV, TZV and QZV 
basis sets. The latter are more efficient for larger sys- 
tems, because more integrals vanish. 



bottleneck in conventional (TD)DFT calculations. RI-J 
replaces the four-index Coulomb integrals by three-index 
and two-index integra ls, w hich considerably lowers the 
algorithmic pre- factor [208j . It is generally safe to use 
with the appropriate auxiliary basis sets. As soon as hy- 
brid functionals arc used, however, the computation of 
the exact exchange becomes rate-determining. 



I. Summary 



G. Diffuse functions 

Adding a [lslpldlf] set of diffuse functions to TZVP 
we obtain the aug-TZVP basis set. The aug-TZVP ex- 
citation energies of all states except the 1 1 A U Rydberg 
state are within 0.01 eV of the reference aug-QZVP re- 
sults and can be considered essentially converged for the 
purposes of present TDDFT. A similar observation can 
be made for the oscillator strengths in Table IIIII 

Going to the even larger quadruple zcta valence (QZV) 
basis sets, the results change only marginally, but the 
computation times increase substantially. In density 
functional theory, these basis sets are mainly used for 
bechmarks and calibration. 



H. Resolution of the identity 

For comparison, we have included results that were ob- 
tained using the resolution of th e identity approximation 
for the Coulomb energy (RI-J) [206l |207( . It is obvious 
that the error introduced by the RI-J approximation is 
much smaller than the basis set error, while the com- 
putation time is reduced by a factor of 5. The RI-J 
approximation is so effective because the computation 
of the Coulomb (Hartree) energy and its respose is the 



To summarize, for larger molecules, SV(P) or similar 
basis sets are often appropriate due to their good cost- 
to-performance ratio. We recommend to check SV(P) 
results by a TZVP calculation whenever possible. Dif- 
fuse functions should be used sparingly for molecules with 
more than about 20 atoms. 



V. PERFORMANCE 

This chapter is devoted to studying and analyzing the 
performance of TDDFT, assuming basis-set convergence. 
We dissect many of the sources of error in typical TDDFT 
calculations. 

To get an overall i mpre ssion, a small survey is given by 
Furche and Ahlrichs [209| |. Typi cal chemical calculations 
are done with the B3LYP 167| functional, and typical 
results arc transition frequencies within 0.4 cV of exper- 
iment, and structural properties of excited states are al- 
most as good as those of ground-state calculations (bond 
lengths to within 1%, dipole moments to within 5%, vi- 
brational frequencies to within 5%). Most importantly, 
this level of accuracy appears sufficient in most cases to 
qualitatively identify the nature of the most intense tran- 
sitions, often debunking cruder models that have been 
used for interpretation for decades. This is proving espe- 
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TABLE IV: Performance of various density functionals for 
the first six singlet excitation energies (in eV) of naphthalene. 
An aug-TZVP basis set and the PBE/TZVP/RI ground-state 
structure was used. The "best" estimates of the true excita- 
tions were from experiment and calculations, as described in 
text. 



Method 


1 1 Ba u 


1 1 B2u 2 1 Ag 


1 ^ig 


2 1 Bs u 


1 M„ 






Pure density functionals 






LSDA 


4.191 


4.026 5.751 


4.940 


5.623 


5.332 


BP86 


4.193 


4.027 5.770 


4.974 


5.627 


5.337 


PBE 


4.193 


4.031 5.753 


4.957 


5.622 


5.141 






Hybrids 








B3LYP 


4.393 


4.282 6.062 


5.422 


5.794 


5.311 


PBEO 


4.474 


4.379 6.205 


5.611 


5.889 


5.603 


"best". 


4.0 


4.5 5.5 


5.5 


5.5 


5.7 



TABLE V: Performance of various wavefunction methods for 
the excitations of Table I. The aug-TZVP basis set and the 
PBE/TZVP/RI ground-state structure was used for a ll ex - 
cept the CASPT2 results, which were ta ken from Ref. |l97| . 
Experimental results are also from Ref. |l97l |. 



Method 
CIS 
CC2 
CASPT2 


1 1 B au 
5.139 
4.376 
4.03 


1 1 B 2u 
4.984 
4.758 
4.56 


2 1 Ag 

7.038 
6.068 
5.39 


1 'Big 

6.251 
5.838 
5.53 


2 x B 3u 
6.770 
6.018 
5.54 


1 U„ 
5.862 
5.736 
5.54 


expt. 


3.97, 4.0 


4.45, 4.7 


5.50, 5.52 


5.28, 5.22 


5.63,5.55 
5.89 




"best" . 


4.0 


4.5 


5.5 


5.5 


5.5 


5.7 



cially useful for the photochemistry of biologically rele- 
vant molecules |69| 



A. Example: Napthalene Results 

As an illustration, compare the performance of var- 
ious density functionals and wavefunction methods for 
the first singlet excited states of naphthalene in Tables 
IIV1 fVl and I VII All calculations were performed using 
the aug-TZVP basis set; the complete active space SCF 
with secon d-ord er perturbation theory (CASPT2) results 
from Ref. |l97| were obtained in a smaller double zeta 
valence basis set with some diffuse augmentation. The 
experimental results correspond to band maxima from 
gas-phase experiments; however, the position of the band 
maximum does not necessarily coincide with the vertical 
excitation energy, especially if the excited state structure 
differs significantly from the ground-state structure. For 
the lower valence states, the CASPT2 results can there- 
fore be expected to be at least as accurate as the experi- 
mental numbers. For higher excited states, the basis set 
used in the CASPT2 calculations appears rather small, 
and the approxi mate second-or der coupled cluster values 
denoted RICC2 [H3, SH H3 might be a better refer- 
ence. Thus our best guess (denoted "best" in the Tables) 
is from experiment for the first 4 transitions, CASPT2 for 
the 5th, and RICC2 for the 6th. 



TABLE VI: Performance of various density functionals and 
correlated wavefunction methods for the oscillator strengths 
of the first three dipole-allowed transitions of naphthalene. 
A aug-TZVP basis set and the PBE/TZVP/RI ground-state 
structure was used for all except the CASPT2 results, which 
were taken from Ref. |l97l |. 



IVlcLIlUU. 


1 1 R, 


1 1 R„ 

1 -D2u 


9 1 R, 


LSDA 


0.0000 


0.0405 


1.1517 


BP86 


0.0000 


0.0411 


1.1552 


PBE 


0.0000 


0.0407 


1.1402 


B3LYP 


0.0000 


0.0539 


1.2413 


PBEO 


0.0000 


0.0574 


1.2719 


LHF/LSDA 


0.0000 


0.0406 


1.2089 


LHF/PBE 


0.0000 


0.0403 


1.2008 


CIS 


0.0002 


0.0743 


1.8908 


CC2 


0.0000 


0.0773 


1.4262 


CASPT2 


0.0004 


0.0496 


1.3365 


expt. 


0.002 


0.102, 0.109 


1.2, 1.3 



We begin with some general observations. 

• The excitation energies predicted by the GGA func- 
tionals BP86 and PBE differ only marginally from 
the LSDA results (an exception being the 1 1 A U Ry- 
dberg state, whose PBE excitation energy is sub- 
stantially lower than those of all other methods). 
Note however that GGA functionals generally im- 
prove over LSDA results for other excited state 
properties such as structures or vibrational frequen- 
cies. 

• Hybrid mixing leads to systematically higher exci- 
tation energies. On average, this is an improvement 
over the GGA results which are systematically too 
low. However, while diffuse excitations benefit from 
hybrid mixing due to a reduction of self-interaction 
error, valence excitation energies are not always im- 
proved, as is obvious for the 1 1 B^ U and 2 1 B% U 
valence states. 

• The 1 1 i?2u state is erroneously predicted below the 
1 1 Bs u state by all density functionals, which is a 
potentially serious problem for applications in pho- 
tochemistry; this is not corrected by hybrid mixing. 

• The configuration-interaction singles (CIS) method 
which uses a Hartree-Fock reference that is compu- 
tationally as expensive as hybrid TDDFT produces 
errors that are substantially larger, especially for 
valence states. The coupled cluster and CASPT2 
methods are far more expensive, and scale pro- 
hibitively as the system size grows. 

The 1 B<2u excitation is polarized along the short 
axis of the naphthalene molecule. In Piatt's nomencla- 
ture of excited states of polycyclic aromatic hydrocar- 
bons (PAHs), 1 1 £?2u corresponds to the 1 L a state. This 
state is of more ionic character than the 1 1 B% U o r 1 L\ > 
state. Parac and Grimme have pointed out [213l | that 
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GGA functionals considerably underestimate the excita- 
tion energy of the 1 L a state in PAHs. This agrees with 
the observation that the 1 1 i?2u excitation of naphtha- 
lene is computed 0.4-0.5 eV too low in energy by LSDA 
and GGA functionals, leading to an incorrect ordering of 
the first two singlet excited states. 



B. Influence of the ground-state potential 

From t he v e ry earliest calculations of transition 
frequencies jl88l Il90l |, it was recognized that the inac- 
curacy of standard density functional approximations 
(LDA, GGA, hybrids) for the ground-state XC potential 
leads to inaccurate KS eigenvalues. Because the approx- 
imate KS potentials have incorrect asymptotic behavior 
(they decay exponentially, instead of as — 1/r, as seen in 
Fig. [J]) , the KS orbital eigenvalues are insufficiently nega- 
tive, the ionization threshold is far too low, and Rydberg 
states are often unbound. 

Given this disastrous behavior, many methods 
have be en deve loped to asymptotically correct 
potentials [2l4 12151 . Any corrections to the ground-state 
potential are dissatisfying, however, as the resulting 
potential is not a functional derivative of an energy func- 
tional. Even mixing one approximation for i> X c( r ) and 
another for / xc has become popular. A more satisfying 
route is to use th e optimized effective potential (OEP) 
method [l59l . Il73| | and include exa ct exchange or other 
self- interaction- free functionals |21fi| . This produces a far 
more accurate KS potential, with the correct asymptotic 
behavior. The chief remaining error is simply the 
correlation contribution to the position of the HOMO, 
i.e, a small shift. All the main features below and just 
above / are retained. 



1. N2, a very small molecule 

A simple system to see the effect of the various ground- 
state potentials is the N2 molecule. In all the cases 
discussed below, a SCF step was carried out using the 
ground-state potential to find the KS levels. These are 
then used as input to Eq. (|47|) with the ALDA XC kernel. 



In Table [VTH the KS energy levels for the LDA 
functional are shown. It is very clear to see that the 
eigenvalues for the higher unoccupied states are positive. 
As mentioned this is due to the LDA potential being 
too shallow and not having the correct asymptotic 
behavior. Comparing the basis-set calculation with the 
basis-set-free calculation, the occupied orbitals are in 
good agreement. However for the unoccupied states 
that are unbounded in LDA, basis sets cannot describe 
these correctly and give a positive energy value which 
can vary greatly from one basis set to another. 



TABLE VII: Orbital energies of the KS energy levels for N2 at 
separation R = 2.0744a.u. Orbitals calculated using the LDA 
potential are shown for two different numerical methods. The 
first is fully numerical basis set free while the other uses the 
Sadlej (52 orbitals) basis set [217| . The OEP method uses the 
EXX (KLI) approximation and is also calculated basis set 
free. 

Energies in (eV) 



Orbital 


LDA 


LDA 


OEP a 




basis set free" 


Sadlej 6 




Occupied orbitals 






la g 


-380.05 


-380.82 


-391.11 


la u 


-380.02 


-380.78 


-391.07 


2a g 


-28.24 


-28.52 


-35.54 


2a u 


-13.44 


-lo.^lU 


90 9Q 




-11.89 


-11.86 


-18.53 


3<7„ 


-10.41 


-10.38 


-17.15 


Unoccupied orbitals 






1tt 9 


-2.21 


-2.23 


-8.44 


4a g 


-0.04 


0.66 


-5.05 


2n u 


> 


1.93 


-4.04 


3a u 


> 


1.35 


-3.54 


is g 


> 




-2.76 


ha g 


> 


3.20 


-2.49 


6a g 


> 




-2.33 


2n g 


> 


3.89 


-2.17 


3tv u 


> 




-2.04 




"From Ref.[218l|. 


'from Ref.[219|. 





TABLE VIII: Comparison of the vertical excitation energies 
for the first twelve excited states of N2 calculated using dif- 
ferent methods for the SCF step. In all cases the KS or- 
bitals from the SCF step are inputed into Casida's equations 
with the ALDA XC kernel. For the LDA calculated with the 
Sadlej basis set, the bare KS transition frequencies are given 
to demonstrate how they are corrected towards their true val- 
ues using Casida's equations. Also given are the mean abso- 
lute errors for each method, errors in backets are calculated 
for the lowest eight transitions only. 









Excitation energy 


(eV) 








State 


Excitation 


BARE KS" ALDA" 


ALDA 6 


LB94" 


OEP rf 


Expt e 








Singlet — ► 


singlet transitions 












— ► iTTg 


9.63 


10.20 


10.27 


9.82 


10.66 


10.27 




Itt „ 


— > lffg 


9.63 


9.63 


9.68 


9.18 


10.09 


9.92 




3<t s 


-> l-Kg 


8.16 


9.04 


9.23 


8.68 


9.76 


9.31 


a Yjg 


3<Jg 


-> 4cr g 






10.48 




12.47 


12.20 




2cr„ 


— > l7Tg 






13.87 




14.32 


13.63 




3cr a 


-» 2tt„ 






11.85 




13.07 


12.90 








Singlet — > 


triplet transitions 








c 3 n u 


2(T„ 


—¥ llTg 


11.21 


10.36 


10.44 


10.06 


11.05 


11.19 


B' 3 £~ 




— > lffg 


9.63 


9.63 


9.68 


9.18 


10.09 


9.67 


W 3 A„ 




— > lWg 


9.63 


8.80 


8.91 


8.32 


9.34 


8.88 


B 3 n B 


3(7g 


— > iTTg 


8.16 


7.50 


7.62 


7.14 


8.12 


8.04 


A 3 E+ 


1tt„ 


— > l7T 9 


9.63 


7.84 


8.07 


7.29 


8.51 


7.74 


E 3 S+ 




— ► 4<Tg 






10.33 


12.32 


11.96 


12.00 


Mean Absolute Error 


(0.61) 


(0.27) 


0.54 


(0.63) 


0.34 





"Using Sadlej basis set. From Ref. [219j|. 
''Basis set free. From Ref. |21SH . 
c From Ref.[22Cj. 

d Using KLI ap proxi mation. From Ref. l21jj|l . 

e Computed in |22ll from the spectroscopic constants of Hubcr 
and Herzberg I222H , 
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In Tabic IVIII1 the bare KS transition frequencies 
between these levels are shown. Note that they are in 
rough agreement with the experimental values and that 
they lie in betw een the singlet-singlet and singlet-triplet 
transitions (223. The ALDA XC kernel f^ DA then 
shifts the KS transitions towards their correct values. 
For the eight lowest transitions LDA does remarkably 
well, the mean absolute error (MAE) being 0.27eV 
for the Sadlcj basis set. For higher transitions it fails 
drastically, the MAE increases to 0.54eV when the next 
four transitions are included. This increase in the MAE 
is attributed to a cancellation of error s tha t lead to good 
frequencies for the lower transitions [218j|. Since LDA 
binds only two unoccupied orbitals, it cannot accurately 
describe transitions to higher orbitals. In basis set 
calculations, the energies of the unbound orbitals which 
have converged will vary wildly and cannot give trusted 
transition frequencies. 

One class of XC functionals that would not have 
this proble m a r e the asymptotica lly corr ected (AC) 
functionals [HI US SH, S2I H^. LB94[227| is one 
such of these and its performance is shown in Table 
IVIIII AC XC potentials tend to be too shallow in the 
core region, so the KS energy levels will be too low 
while the AC piece will force the higher KS states to be 
bound and their energies will cluster below zero. Thus 
it is expected that using AC functionals will consistently 
underestimate the transitions frequencies. 



TABLE IX: Naphthalene: Effect of ground-state potential on 
the excitations of Table HVl A ground-state calculation using 
exact exchange OEP (LHF) is performed and the excitations 
are found using a LDA/PBE kernel respectfully. The result 
is then compared to that found if the LDA/PBE functional 
had been used for both steps. 



Method 


1 1 B$ U 


1 1 B-2u 


2 1 A g 


1 1 Bi a 


2 1 i?3ti 


1 X A U 


LSDA 


4.191 


4.026 


5.751 


4.940 


5.623 


5.332 


LHF/LSDA 


4.317 


4.143 


5.898 


5.097 


5.752 


5.686 


PBE 


4.193 


4.031 


5.753 


4.957 


5.622 


5.141 


LHF/PBE 


4.295 


4.121 


5.876 


5.091 


5.741 


5.693 


"best". 


4.0 


4.5 


5.5 


5.5 


5.5 


5.7 



LHF potential does not change the results significantly. 

The 1 l Bi g and especially the 1 1 A U states are dif- 
fuse, and it is not surprising that their excitation energy 
is considerably underestimated in the LSDA and GGA 
treatment. Using the asymptotically correct LHF poten- 
tial corrects the excitation energy of the 1 A u , which 
is a pure one-particle excitation out of the la u valence 
into the 10a g Rydberg orbital; the latter may be viewed 
3s orbital of the CioHg ion. On the other hand, 
a strong mixture of valence and Rydberg excitations oc- 
curs in 1 1 Bi g . The LHF potential improves the GGA 
results only marginally here, suggesting that more ac- 
curate XC kernels are necessary to properly account for 
valence-Rydberg mixing. 



A much better approach is using the OEP method. 
The KS orbitals found using this method are self- 
interaction free and are usually better approximations 
to the true KS orbitals. OEP will also have the correct 
asymptotic behavior and as we can see in Table IVlTI all 
orbital energies are negative. In Table IVIIII the MAE 
for OEP is 0.34eV, much lower than LDA. Since OEP 
binds all orbitals, it allows many more transitions to be 
calculated. A common OE P fun ctional is exact exchange 
(or the KLI approximation [228j to it) which neglects cor- 
relation effects, but these are generally small contribu- 
tions to the KS orbitals. Using these with ALDA for / xc 
(which does contain correlation) leads to good transition 
frequencies as shown in Table IVIIII Although LDA is 
sometimes closer to the experimental values for the lower 
transitions, the value of OEP lies in it ability to describe 
both the higher and lower transisions. 



2. Napthalene, a small molecule 

Returning to our benchmark case of Naphthalene, us- 
ing more accurate LHF exchange-only potentials from 
Sec. IIIBI together with an LSDA or PBE kernel produces 
excitation energies in between the GGA and the hybrid 
results, except for the 1 1 A u Rydberg state, whose excita- 
tion energy is significantly improved. Whether the LSDA 
kernel or the PBE GGA kernel is used together with an 



C. Analyzing the influence of the XC kernel 

In this section, we discuss the importance of the XC 
kernel in TDDFT calculations. As mentioned earlier, the 
kernels used in practical TDDFT are local or semi-local 
in both space and time. Even hybrids are largely semi- 
local, as they only mix in 20 — 25% exact exchange. 

In realistic calculations, both the ground-state XC po- 
tential and TDDFT XC kernel are approximated. A sim- 
ple way to separate the error in the XC kernel is to look 
at a test case where the exact KS potential is known. 
Figure [5] shows the spectrum of He using the exact KS 
pot entia l, but with the ALDA XC kernel. It does rather 
well (229j (very well, as shall see later in section [VTT when 
we examine atoms in more detail). Very similar results 
are obtained with standard GGA's. 

The errors in such approximate kernels come from the 
locality in space and time. We can test one of these sep- 
arately for the He atom, by studying the exchange limit 
for the XC kernel. For two spin-unpolarizcd electrons, 
f x = — l/2|r — r'|, i.e., it exactly cancels half the Hartree 
term. Most importantly, it is frequency-independent, so 
that there is no memory, i.e., the adiabatic approxima- 
tion is exact. In Fig. [6j we compare ALDAx, i.e., the 
ALDA for just exchange, to the exact exchange result 
for He. Clearly, ALDA makes noticeable errors relative 
to exact exchange, showing that non-locality in space can 
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Continuum 



6S 

24 - 5S- 

4S- 

> 23 - 3P 

«J 3S- 

OJ 
'5b 

S 22 - 



- 




2S 



Exact KS ALDAxc Exact TDDFT 



FIG. 5: T he sp ectrum of Helium calculated using the ALDA 
XC kernel 229] with the exact KS orbitals. 



D. Errors in potential vs. kernel 

In this section, we examine the relative importance of 
the potential and kernel errors. It has long been be- 
lieved that fixing the defects in the potential, especially 
its asymptotic behavior, has been the maj or challenge 
to improving TDDFT results [HI HH, [Hal- We argue 
here that this is overly simplistic, and is due to tests be- 
ing carried out on atoms and small molecules. In large 
molecules, where the interest is in the many low-lying 
transitions, the potential can be sufficiently accurate, 
while the kernel may play a larger role. 

In fact, our analysis of the general failure of TDDFT 
in underestimating the L a transitions in PAH's sheds 
some light on its origin. Using the self-interaction free 
LHF potential does not cure this problem, as is obvious 
from Tab. |IXJ To the best of our knowledge, the cause 
of this shortcoming of TDDFT is not well understood. 
We note, however, that the same incorrect ordering of 
1 L a and 1 Lf ) occurs in the CIS approximation, which is 
self-interaction free. The analysis here shows that this is 
a failure of our approximations to the XC kernel rather 
than to the ground-state potential. 



E. Understanding linear response TDDFT 



Continuun 



6S- 
24 - 5S- 



> 23 



3 22 ■ 



. sp- 
ss- 



FIG. 6: The spectrum of Helium calcu lated using the ALDAx 
kernel and the exact exchange kernel [229T ]. Again the exact 
KS orbitals were used. The importance of non-locality for the 
XC kernel can be seen as the exchange part of ALDA gives 
a noticeable error compared to the exchange part of the true 
functional (the AEXX kernel for He). 



be important. 

Thus the hybrid functionals, by virtue of mixing some 
fraction of exact exchange with GGA, will have only 
slightly different potentials ( mostly in the asymptotic 
region), but noticeably different kernels. 



Several simple methods have evolved for qualitatively 
understanding TDDFT results. The most basic i s the 
single-pole approximation (SPA), which originated |l88j 
in including only one pole of the response function. The 
easiest way to see this here is to truncate Eq. (jT7|) to a 
lxl matrix, yielding an (often excellent) approximation 
to th e change in transition frequency away from its KS 
valuefl93ll230t: 



w- a + 2oj qa K q<jq<J (SPA), 



(58) 



(The original SPA was on the unsymmetric system yield- 
K qaqa , wh ich fo r a spin-saturated system 



"go- 



ing U) 

becomes uj m uj q + 2K qq [2291 ]). This can also provide a 
quick and dirty estimate, since only KS transitions and 
one integral over f xc arc needed. While it allows an es- 
timate of the shift of transitions from their KS energy 
eigenvalue differences, it says nothing about oscillator 
strengths, which are unchanged in SPA from their KS 
values. In fact, a careful analysis of the TDDFT equa- 
tion shows that oscillator strengths are particularly sensi- 
tive to even small off-diagonal matrix elements, whereas 
transition frequencies are less so [231 1. 

A more a dvan ced analysis is the double pole 
approximation [232j (DPA), which applies when two tran- 
sition are strongly coupled to one another, but not 
strongly to the rest of the transitions. Then one can show 
explicitly the very strong effect that off-diagonal elements 
have on oscillator strengths, showing that sometimes an 
entire peak can have almost no contribution. One also 
sees pole-repulsion in the positions of the transitions, a 
phenomenon again missing from SPA. 



TABLE X: Transition frequencies and oscillator strengths 
(O.S) calculated using the double pole approximation (DPA) 
for the lowest transitions in Naphthalene. The PBE 

functional was used with a aug-TZVP basis set on top of a 
PBE/TZVP/RI ground state structure 







KS 


DPA 


Full TDDFT 




LO 


O.S 


LO O.S 


LO O.S 




4.117 


(1.02) 


4.245 (0.001) 


4.191 (0) 


2 Bin 


4.131 


(1.00) 


6.748 (2.02) 


5.633 (1.14) 



The DPA was used recently and very successfully to 
explain X-ray edge spectroscopy results for 3rf-tran sition 
metal solids as one moves across the periodic table |233j . 
These transitions form a perfect test case for DPA, as the 
only difference between them is caused by the spin-orbit 
splitting (several eV) of the 2^> 1 / 2 and 2p 3 / 2 levels. In a 
ground-state KS calculation, this leads to a 2:1 branching 
ratio for the two peaks, based simply on degenearcy, as 
all matrix elements arc identical for the two transitions. 
Experimentally, while this ratio is observed for Fe, large 
deviations occur for other elements. 

These deviations could be seen in full TDDFT calcu- 
lations, and were attributed to strong core-hole correla- 
tions. The SPA, while it nicely accounts for the shifts 
in transition frequencies relative to bare KS transitions, 
but yields only the ideal 2:1 branching ratio. However, 
the DPA model gives a much simpler, and more benign 
interpretation. The sensitivity of oscillator strengths to 
off-diagonal matrix elements means that, even when the 
off-diagonal elements are much smaller than diagonal el- 
ements (of order 1 eV) , they cause rotations in the 2-level 
space, and greatly alter the branching ratio. Thus a KS 
branching ratio occurs even with strong diagonal 'corre- 
lation', so long as off-diagonal XC contributions are truly 
negligible. But even small off-diagonal 'correlation', can 
lead to large deviations from KS branching ratios. 

We can use DPA to understand the lowest 1 B^, U tran- 
sitions in our naphthalene case. In Table IXl we list the 
TDDFT matrix elements for the PBE calculation for the 
two nearly degenerate KS transitions, la u — > 2&3 g and 
26i« — > 262 g , along with their corresponding KS transi- 
tion frequencies. Contour plots of the four orbitals in- 
volved are shown in Fig. \7\ We note first that these 
two KS transitions are essentially degenerate, so that 
there is no way to treat them within SPA. The degen- 
eracy is lifted by the off-diagonal elements, which cause 
the transitions to repel each other, and strongly rotate 
the oscillator strength between the levels, removi ng al - 
most all the oscillator strength from the lower peakj232|. 
The DPA yields almost the correct frequency and oscil- 
lator strength (i.e., none) for the lower transition, but 
the higher one is overestimated, with too much oscil- 
lator strength. This must be due to coupling to other 
higher transitions. In the DPA, in fact the higher transi- 
tion lands right on top of the third transition, so strong 
coupling occurs there too. This example illustrates (i), 
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FIG. 7: The four orbitals involved in the hrst two B^ u 
(contour value ± 0.07 a.u.). The PBE functionals and an 
aug-TZVP basis set were used. 



that solution of the full TDDFT equations is typically 
necessary for large molecules which have many coupled 
transitions, but also (ii), that simple models can aid 
the interpretation of such results. All of which shows 
that, while models developed for well-separated transi- 
tions might provide some insight for specific transitions 
in large molecules, the number and density of transitions 
make such models only semi-quantitative at best. 



VI. ATOMS AS A TEST CASE 

In this section, we look more closely at how well 
TDDFT performs for a few noble gas atoms. As ex- 
plained above, this is far from representative of its be- 
havior for large molecules, but this does allow careful 
study of the electronic spectra without other complica- 
tions. Most importantly, for the He, Be, and Ne atoms, 
we have essentially exact ground-state KS potentials from 
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FIG. 8: Singlet energy level diagram for the helium atom. 
The Rydberg series of transition frequencies clustered below 
the ionization threshold can be seen. The frequencies clus- 
ter together, making it difficult to assess the quality of the 
TDDFT calculated spectra. As discussed in the text, the 
quantum defect is preferable for this purpose. 



Umrigar and coworkers [l56l 1234], This allows us to dis- 
sect the sources of error in TDDFT calculations. 



A. Quantum defect 

In Fig. [5] we show the KS orbital energy level diagram 
of the helium atom. The zero is set at the onset of the 
continuum and is marked with a dotted line. For closed 
shell atoms and for any spherical one-electron potential 
that decays as —1/r at large distances, the bound-state 
transitions form a Rydberg series with frequencies: 



1 



2(n - [inif 



(59) 



TABLE XI: Transition energies (lo) and quantum defects 
(QD) for He atom s-Rydberg series[a.u.]. The ionization en- 
ergy is 0.90372 a.u. 



Transition 


Sing 


;let a 


Triplet a 










LO 


QD 


LO 


QD 


LO 


QD 


Is 




2s 


0.7578 


0.1493 


0.7285 


0.3108 


0.7459 


0.2196 


Is 




3s 


0.8425 


0.1434 


0.8350 


0.3020 


0.8391 


0.2169 


Is 




4s 


0.8701 


0.1417 


0.8672 


0.2994 


0.8688 


0.2149 


Is 




5s 


0.8825 


0.1409 


0.8883 


0.2984 


0.8818 


0.2146 


Is 




6s 


0.8892 


0.1405 


0.8926 


0.2978 


0.8888 


0.2144 


Is 




7s 


0.8931 


0.1403 


0.8926 


0.2975 


0.8929 


0.2143 



"Accurate non-relativistic calculations from Ref. |19q |. 
b The differences betwee n the KS eigenvalues obtained with the 
exact potential from Ref. I156H . 
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FIG. 9: The exact s KS quantum defect and the exact singlet 
and triplet quantum defects of He and their parabolic fits. 
The quantum defect may clearly be described as a smooth 
function of energy, in this linear fit. Thus knowing the 

quantum defect for a few transitions allows us to find it for 
all transitions, and hence their frequencies. 



where / is the ionization potential, and (i n i is called 
the quantum def ect. Quantum d efect theory was de- 
veloped by Ham [235| | and Seaton |236| before even the 
Hohcnberg-Kohn theorem[l|. 

The great value of the quantum defect is its ability to 
capture all the information about the entire Rydberg se- 
ries of transitions in a single slowly-varying function, the 
quantum defect as a function of energy, (ii(E = lo — I), 
which can often be fit by a straight line or parabola. 
In Table IXIl we report extremely accurate results from 
wavefunction calculations for the helium atom. We show 
singlet and triplet values that have been obtained by 
Drake [198| | . We also give results fro m th e exact ground- 
state KS potential shown in Fig [3] |156| . For each col- 
umn, on the left are the transition frequencies, while on [5] 



the right are the corresponding quantum defects. Note 
how small the differences between transitions become as 
one climbs up the ladder, and yet the quantum defect 
remains finite and converges to a definite value. 

All the information of the levels of Fig. [8] and of Ta- 
ble |XI] is contained in Fig. This clearly illustrates that 
the quantum defect is a smooth function of energy, and is 
well approximated (in these cases) as a straight line. The 
quantum defect is thus an extremely compact and sensi- 
tive test of approximations to transition frequencies. Any 
approximate ground-state KS potential suggested for use 
in TDDFT should have its quantum defect compared 
with the exact KS quantum defect, while any approxi- 
mate XC-kernel should produce accurate corrections to 
the ground-state KS quantum defect, on the scale of Fig. 
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FIG. 10: The Be p quantum defect of LB94, exact exchange 
(OEP), and KS, and their best fits. While both functionals 
give the correct asymptotic behavior of the KS potential, by 
calculating the quantum defect, we can learn more about their 
performance. 



To demonstrate the power of this analysis, we test 
two common approximations to the ground-state poten- 
tial, both of which produce asympt otica lly correct po- 
tentials. T hese are exact exchange [237| (see Sec III B|) 
and LB94 |227| . Exact exchange calculations are more 
demanding than traditional DFT calculations, but are 
becomi ng popula r because of the high quality of the po- 
tential |238l . [239[. On the other hand, LB94 provides 
an asymptotically correc t po t ential at little extra cost 
beyond traditional DFT [Hi HH, [HI- In Fig. [TO] we 
show the p Be quantum defect obtained with LB94, OEP, 
and exact KS potentials. Fig. [TOI immediately shows the 
high quality of the exact exchange potential. The quan- 
tum defect curve is almost identical to the exact one, 
apart from being offset by about 0.1 . On the other 
hand, t he quantu m defect of LB94 was poor for all cases 
studied [24ll . |242| . This shows that just having a poten- 
tial that is asymptotically correct is not enough to get a 
good quantum defect. 



B. Testing TDDFT 

To see how well TDDFT really does, we plot quantum 
defects for atoms. We take the He atom as our proto- 
type, as usual in this section. In Fig. 111! we plot first 
the KS quantum defect and the exact singlet and triplet 
lines, as before in Fig. [9] Then we consider the Hartrec 
approximation. This is equivalent to setting the XC ker- 
nel to zero. This changes the postion of the singlet curve, 
but leaves the triplet unchanged from its KS value, be- 
cause the direct term includes no spin-flipping. It defi- 
nitely improves over the KS for the singlet. Lastly, we 
include ALDA XC corrections. Only if these significantly 
improve over the Hartrec curves can we say TDDFT is 



o 

CD 

a5 0.3 







O 


ALDA 


0.4 




A 


Hartree 








KS 






• 


Exact 



CO 
=! 

o 



0.2 



0.1 



0.0 



-• - 

— Q- 



-e — o-e — 



- — e-e- 



-A — frA- 



-0.15 -0.10 -0.05 

Energy [a.u.] 



0.00 



FIG. 11: The corrections due to using the Hartree or ALDA 
kernel on the exact KS s quantum defect of He. Using the 
Hartree kernel only effects the singlet values, shifting them too 
low. If a good XC kernel is then used, it should move both the 
triplet and singlet qua ntum defects from the Hartree kernel 
towards the exact ones [242l |. In this case, ALDA does a good 
job and is performing well. 



really working here. Clearly it does, reducing the Hartree 
error enormously. 

These results are also typical of He P transitions, and 
Be S and P transitions. For reasons as yet unclear the 
s — > d transitions fail badly for both these systems [24ll 
l24l. 



C. Saving standard functionals 

We have a problem with the incorrect long-range be- 
havior of the potential from standard density functionals 
only when Rydbcrg excitations are needed. But it would 
be unsatisfactory to have to perform a completely differ- 
ent type of calculation, eg OEP, in order to include such 
excitations when desired, especially if the cost of that 
calculation is significantly greater. 

However, it is possible, with some thought and care, 
and using quantum defect theory, to extract the Ryd- 
berg series from the short-ranged LDA potential. To see 
this, consider Fig. [12l which shows both the bare KS 
response and the TDDFT corrected response of the He 
atom. The <5-function absorptions at the discrete transi- 
tions have been replaced by straightlines, whose height 
represents the oscillator strength of the abs orpti on, mul- 
tiplied by the appropriate density of states 2471 ] . In the 
top panel, just the KS transitions are shown, for both 
the KS potential and the LDA potential of Fig [3] from 
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FIG. 12: He atom: The top panel shows the bare exact KS 
and LDA spectra, and the lower panel shows the TD DFT 
corrected spectra, LDA/ALDA results a re fr om [244 [ but 
unshifted; the exact calculations are from [245| . multiplied by 
the density of st ates factor (see text), and the experimental 
results are from [2491 



section Hi Bl The exact curve has a Rydberg series con- 
verging to 0.903, the exact ionization threshold for He. 
The LDA curve, on the other hand, has a threshold at 
just below 0.6. But clearly its optical absorption mimics 
that of the exact system, even in the Rydberg series re- 
gion, and is accurate to about 20%. The TDDFT ALDA 
corrections are small, and overcorrect the bare LDA re- 
sults, but clearly are consistent with our observations for 
the bare spectra. 

Why is this the case? Is this a coincidence? Returning 
to Fig. [3] of the introduction, we notice that the LDA 
(or GGA) potential runs almost exactly parallel to the 
true potential for r < 2, i.e., where the density is. Thus 
the scattering orbitals of the LDA potential, with transi- 
tion energies between 0.6 and 0.9, almost exactly match 
the Rydberg orbitals of the exact KS potential with the 
same energy. When carefully defined, i.e., phase space 
factors for the continuum relative to bound states, the 
oscillator strength is about the same. This is no coinci- 
dence but, due to the lack of derivative discontinuity of 
LDA, its potential differs from the exact one by roughly 
a constant. 

The 'fruitfly' of TDDFT benchmarks is the tt -> tt* 
transition in benzene. This occurs at about 5 eV in a 
ground-state LDA c alcula tion, and ALDA shifts it cor- 
rectly to about 7 eV [230j |. Unfortunately, this is in the 
LDA continuum, which starts at about 6.5 eV! This is an 
example of the same general phenomenon, where LDA 
has pushed some oscillator strength into the continuum, 
but its overall contribution remains about right. 

We can go one step further, and even deduce the en- 
ergies of individual transitions. While the existence of a 
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FIG. 13: He atom: solution of Eq. (|60[l for jjb as a function of 
r; The n = 20 orbital was used for the exact case, and the 
scattering orbital or energy E — I + e%P A was used for the 
LDA. 



quantum defect requires a long-ranged potential, its value 
is determined by the phase-shift caused by the deviation 
from — 1/r in the inte rior of the atom. The quantum 
defect extractor (QDE)[248|, is a formula for extracting 
the effective quantum defect from a scattering orbital of 
a short-ranged KS potential, such as that of LDA. The 
QDE is: 



dliK, 
dr 



1 



n* 1 U(-n*;2;2r/n*) 
T rU(l - n*;2;2r/n*) 



(60) 



Here k = \/2\E\ is written as k — (n*)^ 1 , with n* — 
(n — /!„), where n numbers the bound state, and fj. n is 
the quan tum defect; U is the confluent hypergeometric 
function [249j | . If the extractor is applied to an orbital of a 
long-ranged potential, it rapidly approaches its quantum 
defect. 

In Fig [13l we plot the results of the QDE for the He 
atom, applied to both the exact KS potential and the 
LDA potential. The LDA potential runs almost paral- 
lel to the exact one in the region 1 < r < 2 (where 
fioo can already be extracted accurately), and orbitals 
corresponding to the same frequency (exact and LDA) 
are there fore very close in that region. In the spirit of 
Ref. [250j . we compare the exact energy-normalized 20s 
orbital (which is essentially identical to the zero-energy 
state in the region < r < 6) and the LDA orbital of 
energy I + e\^ A = 0.904 - 0.571 = 0.333. Notice how 
good the LDA orbital is in the region 1 < r < 2. We 
show in FigfT^] the solution of Eq. (|6"0"|) when this scat- 
tering LDA orbital is employed. Clearly, the plateau of 
the LDA curve in the 1 < r < 2 region is an accurate 
estimate of the quantum defect. The value of fi on this 
plateau is 0.205, an underestimation of less than 4% with 
respect to the exact value. 

Thus, given the ionization potential of the system, 
LDA gives a very accurate prediction of the asymptotic 
quantum defect. The ionization potential is needed to 
choose the appropriate LDA scattering orbital, but the 
results are not terribly sensitive to it. We repeated the 
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same procedure with the LDA ionization potential (de- 
fined as £ LDA (He)-£: LDA (He + )=0.974) instead of the 
exact one, and found £t^ A = 0.216, overestimating the 
exact /loo by just 1%. 



D. Electron scattering 

Lastly in this section, we mention recent progress in de- 
veloping a theory for low-energy electron scattering from 
molecules. This was one of the original motivations for 
developing TDDFT. One approach would be to evolve a 
wavcpacket using the TDKS equ ations, bu t a more di- 
rect approach has been developed [25ll \25^ , in terms of 
the response function \ °f the N + 1 electron system 
(assuming it is bound). 

This uses similar technology to the discrete transition 
case. Initial results for the simplest case, electron scat- 
tering from He + , suggest a level of accuracy compara- 
ble to bound-bound transitions, at least for low energies 
(the most difficult cas e for traditional methods, due to 
bound-free correlation [253j h TDDFT, using the exact 
ground-state potential and ALD A, pro duces more accu- 
rate answers than static exchange [2 5 1| . a tra ditional lo w 
cost method that is used for larger molecules [254 [255j . 

However, that TDDFT method fails when applied to 
electron scattering from Hydrogen, the true prototype, 
as the approximate solution of the TDDFT equations 
(very similar to the single pole approximation of Sec IV Ep 
fails, due to stronger correlations. To overcome this, a 
much simpler met hod h as been developed, that uses an 
old scattering trick |256| to deduce scattering phase shifts 
from bound-state energies when the system is placed in a 
box, yielding excellent results for a very demanding case. 



VII. BEYOND STANDARD FUNCTION ALS 

We have surveyed and illustrated some of the many 
present successful applications of TDDFT in the previous 
section. In these applications, standard approximations 
(local, gradient-corrected, and hybrid, see sec. Ill B[) arc 
used both for the ground-state calculation and the exci- 
tations, via the adiabatic approximation (sec. IIII D[) . In 
this section, we survey several important areas in which 
this approach has been found to fail, and what might be 
done about it. 

The errors are due to locality in both space and time, 
and these are intimately related. In fact, all me mory 
effects, i.e., dependence on the history of the density jl84| . 
implying a frequency-dependence in the XC k erne l, can 
be subsumed into an initial-state dependence [HO], but 
probably not vice-versa. Several groups are attempting 
to build such e ffect s into new approximate f uncti onals 

[257L m, m sn hed, eh hm hh hm hh, but 



and foun d to be missing in the usual linear response 
treatment [l90| |. The failure of the local approximation 
in space is seen when TDDFT is applied to extended 
systems, e.g., polymers or solids. Local approxima- 
tions yield short-ranged XC kernels, which become ir- 
relevant compared to Hartrce contributions in the long- 
wavelength limit. The Coulomb repulsion between elec- 
trons generally requires long-ranged (ie 1/r) exchange ef- 
fects when long-wavelength response is being calculated. 

Thus several approaches have been developed and ap- 
plied in places where the standard formulation has failed. 
These approaches fall into two distinct categories. On 
the one hand, where approximations that are local in 
the density fail, approximations that are local (or semi- 
local) in the current- density might work. In fact, for 
TDDFT, the gradient expansion, producing the leading 
corrections to A LDA, only works if the current is the ba- 
sic variable |267l |. Using the gradient expansi on itself i s 
called the Vignale-Kohn (VK) approximation [268l |269| |. 
and it has been tried on a variety of problems. 

The alternative approach is to construct orbital- 
dependent app r oxim ations with explicit frequency- 
dependence [270, l27ll |. This can work well for specific 
cases, but it is then hard to see how to construct general 
density functional approximations from these examples. 
More importantly, solution of the OEP equations is typ- 
ically far more expensive than the simple KS equations, 
making OEP impractical for large molecules. 



Double excitations 



none have shown universal applicability yet. 

The failure of the adiabatic approximation is most no- 
ticeable when higher-order excitations are considered, 



As first pointed out by Casida[190j, double excita- 
tions appear to be missing from TDDFT linear resp onse , 
with in any adiabatic approximation. Experience [2 72l 
|273| ] shows that, like in naphthalene, sometimes adia- 
batic TDDFT will produce a single excitation in about 
the right region, in place of two lines, where a double has 
mixed strongly with a single. 

In fact, when a double excitation lies close to a single 
excitation, elementary quantum mechanics s hows that 
/ xc must have a strong frequency dependence [27C| . Re- 
cently, po st-ad i abat i c TD DFT methodologies have been 
developed [2 70l 1274 l275l | for including a double excita- 
tion when it is close to an optically- active sin gle excita- 
tion, and works well for small dienes [270l l276j . It might 
be hoped that, by going beyond linear response, non- 
trivial double excitations would be naturally included in, 
e.g., TDLDA, but it has recently been proven that, in 
the higher-order response in TDLDA, the double excita- 
tions occur simply at the sum of single-excitations [277j | . 
Thus we do not currently know how best to approximate 
these excitations. This problem is particularly severe for 
quantum wells, where the external potential is parabolic, 
leading to mu ltiple near degeneracies between levels of 
excitation 274 1 . 

Returning to our naphthalene example, based on a HF 
reference, the 2 1 A g state has, according to the RICC2 
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results, a considerable admixture of double excitations. 
This is consistent with the fact that the CIS method 
yields an excitation energy that is too high by 1.5 eV 
compared to experiment. The TDDFT results are much 
closer, yet too high by several tenths of eV. 

B. Polymers 



experiment [284 12851 . 12861. An alternative is to study 
the many-body problem [2871 ]. and ask which expressions 
must the XC kernel i nclude in order to yield an accurate 
absorption spectrum [288l . [289| . However, the presently 
available scheme s req uire an expensive GW calculation 
in th e firs t place [290j ] . A recent review can be found in 
Ref. HEtt 
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FIG. 14: ALDA and VK static axial polarizability of poly- 
acetylene com pared with RHF and MP2 results from Refs. 
[2781 . |279| . |280| |. ALDA severely overestimates the polarizabil- 
ity compared to the accurate MP2 calculation. Hartree-Fock 
is also incorrect. However using the VK functional gives al- 
most exact agreement, at least in this case. 



D. Charge transfer 

As is usually the case whenever a method is shown to 
work well, it starts being applied to many cases, and spe- 
cific failures appear. Charge transfer excitations are of 
great importance in photochemistry, especially of biolog- 
ical systems, but many workers have now found abysmal 
results with TDDFT for these cases. 

This can be understood from the fact that TDDFT 
is a linear response theory. When an excitation moves 
charge from one area in a molecule to another, both ends 
will relax. In fact, charge transfer between molecules can 
be well-approximated by ground-state density functional 
calculations of the total energies of the species involved. 
But TDDFT must deduce the correct transitions by in- 
finitesimal perturbations around the ground-state, with- 
out an relaxation. Thus it seems a poor problem to tackle 
with linear response. Many researchers are studying this 
problem, to understand it and find practical solutions 
around it [H, [293, SH Hi, [H| . 



An early triumph of the VK functional was the static 
polarizabilitics of long-chain conjugated polymers. These 
polarizabilites are greatly underestimated by LDA or 
GGA, with the error growing rapidly with the number 
of units [28l| . On the other hand, HF does rather well, 
and does not overpolarize. The VK correction to LDA 
yields excellent results in many (but not all) cases, show- 
ing that a current-dependent functional can correct the 
over-polarization problem. Naturally, orbi tal-dependent 
functionals also account for this effect [2821 ] . but at much 
higher computational cost. 

C. Solids 

Again, in trying to use TDDFT to calculate the opti- 
cal response of insulators, local approximations has been 
shown to fail ba dly. M ost noticeably, they do not describe 
excitonic effects [283j, or the exciton spectrum within the 
band gap. On top of this, the gap is usually much smaller 
than experiment, because adiabatic appoximations can- 
not change the gap size from its KS value. 

One approach is using the VK approximation in TD- 
CDFT. This has proven rather successful, although a sin- 
gle empirical factor was needed to get agreement with 



VIII. OTHER TOPICS 

In this chapter, we discuss several topics of specialized 
interest, where TDDFT is being applied and developed 
in ways other than simple extraction of excitations from 
linear response. 

In the first of these, we show how TDDFT can be used 
to construct entirely new approximations to the ground- 
state XC energy. This method is particularly useful for 
capturing the long-range fluctuations that produce dis- 
persion forces between molecules, which are notoriously 
absent from most ground-state approximations. 

In the second, we briefly survey strong field applica- 
tions, in which TDDFT is being used to model atoms 
and molecules in strong laser fields. We find that it works 
well and easily for some properties, but less so for others. 

In the last, we discuss the more recent, hot area of 
molecular electronics. Here, many workers are using 
ground-state DFT to calculate transport characteristics, 
but a more careful formulation can be done only within 
(and beyond) TDDFT. We review recent progress toward 
a more rigorous formulation of this problem. 
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FIG. 15: Binding energy for the Helium dime r int eracting 
via Van der Waals (VdW) forces, from Ref. [296l ]. Using 
the fluctuation-dissipation theorem (FDT), new XC energy 
functionals may be constructed using any ground-state func- 
tional. The c urves from the standard ground-state functionals 
BP86[!,[297| and PBE0 are given, as well as the FDT- curves 
with these as input. Clearly the FDT is needed to accurately 
describe VdW interaction. 



A. Ground-state XC energy 



TDDFT offers a method[296j, \2M, |299|, l30fj, 1301| to 
find more sophisticated ground-state approximate en- 
ergy functional using the frequency-dependent response 
function. Below we introduce the basic formula and 
discuss some of the exciting systems this method is 
being used to study. 

This procedure uses the adiabatic connection 
fluctuation-dissipation formula: 



E xc [n ] = lJ o d\ J d\ J d 



(61) 



where the pair density is 



Pl(r, r') 




Xlaa' [no]{rr';iw) 



-n (r)6 



(3), 



and the coupling-constant A is defined to multiply the 
electron-electron repulsion in the Hamiltonian, but 
the external p otential is adjusted to keep the density 
fixed[30i, Ho3. xlaa' is given by Eq (42]) with the XC 
kernel /Z X c<t<t' ■ Any approximation to the XC kernel 
yields a sophisticated XC energy E xc [n]. 

It is interesting that if we set XC effects to zero in 
conventional DFT, we end up with the highly inaccurate 



Hartree method of 1928. However when calculating the 
linear response, if the XC kernel is zero (i.e. within the 
random phase approximation) , the XC energy calculated 
using Eq. (|61[) still gives useful results. 



Computationally this procedure is far more demand- 
ing than conventional DFT, but as the above example has 
shown, even poor approximations to the XC kernel can 
still lead to good results. Using this method to find the 
XC energy has the ability to capture effects such as dy- 
namical correlation or Van der Waal interactions, which 
are missing from conventional ground-state DFT approx- 
imations, and are thought to be important in biological 
systems. 

In particular, the coefficient in the decay of the energy 
between two such pieces (Cfe in E — > —Ca/R e , where 
R is their separation) can be accurately (within about 
20%) evaluated using a local a pproximat i on to th e 
frequency-dependent polarizabilitv [299T |304| . 13051 . l306l |. 
In Fig. 1151 the binding energy curve for two Helium 
atoms interacting via Van der Waals is shown. Using 
the fluctuation-dissipation formula, Eq. (|6ip . and the 
PBEO XC kernel clearly gives more accurate results 
than semi- local functionals. Recently, the frequency 
integral in Eq. (|6ip has been performed explicitly but 
approxim ately , yie l ding an explic i t non-local density 
functional [293, [2M HqI HhI [3(1, [H, EH EH E3 
applicable at all separations. TDDFT response functions 
have also been used in the framework of symmetry- 
adapted perturbation theory to generate acc urat e 
binding energy curves of Van der Waals molecules |313| . 



One can go t he other w ay, and try using Eq. (|6ip for 
all bond lengths [HI EH- In fact > E Q © provides a KS 
density functional that allo ws b ond-breaking without 
artifical symmetry breaking [300l |. In the paradigm 
case of the H2 molecule, the binding energy curve has 
no Coulson-Fischer point, and the dissociation occurs 
correctly to two isolated H atoms. Unfortunately, 
simple approximations, while yielding correct results 
near equilibrium and at infinity, produce an unphysical 
repulsion a t lar ge but finite separations. This can be 
traced back [300j | to the lack of double excitations in any 
adiabatic f xc . Study of the convergence of E KC with 
basis sets has also led t o an obvious flaw in the ALDA 
kernel at short distances 12961. 



Further w ork is needed to find accurate XC kernels. 
One method [298| to test these is by examining the uni- 
form electron gas as the frequency dependend susceptibil- 
ity can be found easily and uses the well known Lindhard 
function. Hence different approximate XC kernels may 
be tested and their results compared to highly accurate 
Monte-carlo simulations. 
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B. Strong fields 

Next we turn our attention to the non-perturbative 
regime. Due to advances in laser technology over 
the past decade, many experiments are now possible 
in regimes where the laser field is stronger than the 
nuclear attraction [ll|. The time-dependent field cannot 
be treated perturbatively, and even solving the time- 
dependent Schrodinger equation in three dimensions 
for the evolution of two interacting electrons is b arely 
feasible with present-day computer technologv [316| . 

For more electrons in a time-dependent field, wavefunc- 
tion methods are prohibitive, and in the regime of (not 
too high) laser intensities, where the electron-electron in- 
teraction is still of importance, TDDFT i s essentially the 
only prac tical scheme available [3171 13181 . 13191 l320l . l32ll 
|322| . |323| . There are a whole host of phenomena that 
TDDFT might be able to predict: high harmonic gen- 
eration, multi-photon ionization, above-threshold ioniza- 
tion, above-threshold dissociation, etc., but only if accu- 
rate approximations are available. 

With the recent advent of atto-second laser pulses, the 
electronic time-scale has become accessible. Theoretical 
tools to analyze the dynamics of excitation processes 
on the attosecond time scale will become more and 
more important. An example of such a tool is the 
time - depen dent electron localization function (TDELF) 
[325l l326l |. This quantity allows the time-resolved 
observation of the formation, modulation, and breaking 
of chemical bonds, thus providing a visual understanding 
of th e dyn amics of excited electrons (for an example see 
Ref. [Ml). The natural way of calculating the TDELF 
is from the TDKS orbitals. 

High harmonic generation (HHG) is the production 
from medium intensity lasers of very many harmon- 
ics (sometimes hundreds) of the input intensity. Here 
TDDFT calcul ations have bee n rather su ccesful for atoms 
[327l 13281 . l320j | and molecules |330l . 13311 ]. Recent experi- 
ments have used the HHG resp onse of molecules to de- 
termine their vibrational modes [332(1 . Calculations hav e 
been performed using traditional scattering theory |333l |. 
If this method grows to be a new spectroscopy, perhaps 
the electron scattering theory of Sec I VI Dl will be used to 
treat large molecules. 

Multi-photon ionization occurs when an atom or 
molecule loses more than one electron in an intense 
field. About a decade ago, this was discovered to be 
a non-sequential process, i.e., the probability of dou- 
ble ionization can be much greater than the product of 
two independent ionization events, leading to a 'knee' 
in the double ionizat ion probability as a function of 
intensity [HI [HI [Hi|. TDDFT calculations have so far 
been unable to accurately reproduce this knee, and it has 
recently been shown that a correlation-induced deriva- 
tive disc ontin uity is needed in the time-dependent KS 
potential jl85| . 



Above-threshold ionization (ATI) refers to the proba- 
bility of ionization when the laser frequency is less than 
the ionizatio n po t enti al, i.e., it does not occur in lin- 
ear response j337l l338l |. Again, this is not well given 
by TDDFT calculations, but both this and MPI require 
knowledge of the correlated wavefunction, which is not 
directly available in a KS calculation. 

Since the ionization threshold plays a crucial role in 
most strong field phenomena, Koopmans theorem relat- 
ing the energy level of the KS HOMO to the ionization 
energy must be well satisfied. This suggest s th e use 
of self-interaction free methods such as OEP [l59l Il73| 
or LDA-SIC rather than the usual DFT approximations 
(LDA,GGA,etc), with their poor potentials (see Fig. [3] 
in Sec. [HBl) . 

The field of quantum control has mainly concentrated 
on the motion of the nuclear wave packet on a given set 
of precalculated potential energy surfaces, the ultimate 
goal being the femto-second control of chemical reactions 
|339l |. With the advent of atto-second pulses, control of 
electronic dynamics has come within reach. A marriage 
of optimal-control theory with TDDFT appear s to be th e 
ideal theoreti cal tool to tack le these problems [34Ct [34lj . 
Recent workHH, [HI [344[ has shown the ability of 
TDDFT to predict the coherent control of quantum wells 
using Terahertz lasers. However they remains many dif- 
ficulties and challen ges, including the coupling between 
nuclei and electrons [345l 13461 l34l| , in order to develop a 
general purpose theory. 



C. Transport 

There is enormous interest in transport through 

single molecules as a key component in future 

nanotechnology |348| . Present formulations use ground- 
state density functionals to describ e the stationary non- 
equilibrium current-carrying state [34^ |. But several 
recent s ugge stions consider this as a time-dependent 
problem[35a, HHL HI [H|, and use TD(C)DFT for a 
full description of the situation. Only time will tell if 
TDDFT is really needed for an accurate description of 
these devices. 

Imagine the setup shown in Fig. nowhere a conducting 
molecule is sandwiched between two contacts which are 
connected to semi-infinity leads. The Landauer formula 
for the current is 



7= I r dE T(E)[f L (E) - f R (E)] 
* J -00 



(62) 



where T(E) is the transmission probability for a given 
energy and /l/r{E) is the Fermi distribution function 
for the left/right lead. The transmission probability can 
be written using the non-equilibrium Green's functions 
(NEGF) of the system. Ground-state DFT is used 
to find the KS orbitals and energies of the extended 
molecule and used to find the self-energies of the leads. 
These are then fed into the NEGF method, which will 
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continuity and self interacting errors (SIE) in the 
approximation s to t he XC functional may be the source 
of the problem 354]. In Hartree-Fock (HF) calculations 
(and also in OPM calculations [282] with EXX, exact 
exchange), which have no SIE, the conductances come 
out a lot lower in most regions 356]. Also calculations 



FIG. 16: Schematic representation of a benzene-l,4-di-thiol 
molecule between two gold contacts. The molecule plus gold 
pyramids (55 atoms each) constitute the extended molecule as 
used in the DFT calculations for the Landauer approach. 



determine T{E) and hence the current. 

The NEGF scheme has had a number of successes, 
most notably for atomic-scale point contacts and metal- 
lic wires. Generally it does well for systems where the 
conductance it high. However, it was found that for 
molecular wires, the conductance is overestimated by 
1 — 3 orders of magnitude. Various explanations for this 
and the problems with DFT combined with NEGF in 
general have been suggested. 

Firstly, the use of the KS orbitals and energy levels 
has no theoretical basis. The KS orbitals are those 
orbitals for the non-interacting problem that reproduce 
the correct ground-state density. They should not be 
thought of as the true single-particle excitations of the 
true system. However as we have seen they often repro- 
duce these excitations qualitatively, so it is not clear to 
what extent this problem affects the conductance. 

The geometry of the molecules was also suggested as a 
source of error. DFT first relaxes the molecule to find its 
geometry, whereas in the experiments the molecule may 
be subject to various stresses that could rotate parts of 
it and/or squash parts together. However calcul ation s 
have shown that the geometry corrections are small (3541 ] . 

The approximation that the non-equilibrium XC 
functional is the same as for the static case has been 
suggested as a major source of error. In fact neither the 
HK theorem nor the RG theorem are strictly valid for 
current-carrying systems in homogeneous electric fields. 
A dynamical correction to the LDA functional for the 
static case has been derived using the Vignale-Kohn 
functional TDCDFT but were found to yield only small 
corrections to ALDAl355l. 



In a similar vein, the lack of the derivative dis- 



have been done using a simple model 3571 ] with a KS 
potential with a derivative discontinuity. The I-V curves 
for this system are significantly different from those 
predicted by LDA. This problem is most severe when 
the molecule is not strongly coupled to the leads, but 
goes away when it is covalently bonded. Recent OEP 
calculations of the tran smission along a H-atom chain 
verify these features [356j . 

Despite these problems, quantitative results can be 
found for molecular devices. By looking at what bias a 
KS energy level gets moved between the two chemical 
potentials of the leads ( and hence by Eq. (|6"2")) there 
should be a conductance peak), o ne c an qualitatively 
predict postitions of these peaks |358| . although the 
magnitude of the conductance may be incorrect by 
orders of magnitude. 

Since transport is a non-equilibrium process, we should 
expect that using static DFT will not be able to accu- 
rately predict all the features. Recently a number of 
methods have been sugg ested to use TDDFT to calculate 
transport. In Rcf. [359], the authors present a practical 
scheme using TDDFT to calculate current. The basic 
idea is to 'pump' the system into a non-equilibrium ini- 
tial state by some external bias and then allow the KS 
orbitals to evolve in time via the TDKS equations. The 
RG theorem then allows one to extract the longitudinal 
current using the continuity equation. Using transparent 
boundary conditions in the leads (these solve problems 
with propagating KS in the semi-infinite leads) and us- 
ing an iterative procedure to get the correct initial state, 
they are able to find the steady state current. 

An alternative formulation uses pe riod ic bound- 
ary conditions and includes dissipation 360] . In the 
Landauer-Biittiker formulism, dissipation effects due 
to electron-electron interaction and electron-phonon 
interaction are neglected as the molecule is smaller 
than the scattering length. However, there would be 
scattering in the leads. Imagine a molecule in the ring 
geometry, with a spatially constant electric field. Via a 
gauge transformation, this can be replaced by a constant 
time-dependent magnetic field through the center of the 
ring. If there is no dissipation, the electrons would keep 
accelerating indefinitely and never reach a steady state. 

In the classical Boltzmann equation for transport, scat- 
tering is included via a dissipation term using r, the av- 
erage collision time. A master equation approach is ba- 
sically a generalization of the Boltzmann equation to a 
fully quantum mechanical system. The master equation 
is based on the Louville equation in quantum mechanics 



27 





t 




J 



FIG. 17: Ring geometry for gauge transformation of electric 
fields. 



and for a quantum mechanical density coupled to a heat 
bath (or reservoir), it is written as 

j t p(t) = -i[H,p(t)} + Cip(t)} (63) 

where C is a superoperator acting on the density whose 
elements are calculated using Fermi's Golden rule with 
V e i-ph in a certain approximation (weak coupling and 
instantaneous processes). A KS master equation [353| 
can be set up, modifying C for single particle reduced 
density matrices so that it will give the correct steady 
state. The TDKS equations are then used to prop- 
agate forward in time until the correct steady state 
density is found. The current in then extracted from 
this. Recent calculations have shown it can give cor- 
rect behaviour, such as hysteresis in I-V curves. [352ll36l| 



IX. SUMMARY 

We hope we have conveyed some of the spirit and ex- 
citement of TDDFT in this non-comprehensive review. 
We have explained what TDDFT is, and where it comes 
from. We have shown that it is being used, and often 



works well, for many molecular excitations. Its useful- 
ness lies neither in high accuracy nor reliability, but in 
its qualitative ability to yield roughly correct absorption 
spectra for molecules of perhaps several hundred atoms. 
Thus we emphasize that, usually, there are many exci- 
tations of the same symmetry, all coupled together, and 
that these are the circumstances under which the theory 
should be tested. For many molecular systems, TDDFT 
is now a routine tool that produces useful accuracy with 
reasonable confidence. 

That said, we have discussed some of the areas where 
TDDFT in its current incarnation is not working, such 
as double excitations, charge transfer, and extended sys- 
tems. But there has been significant progress in two 
out of three of these, both in understanding the origin 
of the problem, and finding alternative approaches that 
may ultimately yield a practical solution. We also stud- 
ied how well TDDFT works for a few cases where the 
exact ground-state solution is known, describing the ac- 
curacy of different functionals. We also surveyed some 
applications beyond simple linear response for optical ab- 
sorption, such as ground-state functionals from the adi- 
abatic connection, strong fields, and transport. In each 
of these areas, more development work seems needed be- 
fore TDDFT calculations can become a routine tool with 
useful accuracy. 

Many wonder how long DFT's preemminence in elec- 
tronic structure can last. For sure, Kohn-Sham DFT is a 
poor player that struts and frets his hour upon the stage 
of electronic structure, and then is heard no more. After 
all, its predecessor, Thomas-Fermi theory, is now obso- 
lete, being too inaccurate for modern needs. Many al- 
ternatives for electronic excitations, such as GW, are be- 
coming computationally feasible for interesting systems. 
But we believe DFT, and TDDFT, should dominate for 
a few decades yet. 

We thank Michael Vitarclli for early work and Dr. 
Meta van Faassen and Dr. Max Koentopp for pro- 
viding figures. K.B. gratefully acknowledges support 
of the US Department of Energy, under grant number 
DE-FG02-01ER45928, and the NSF, under grant CHE- 
0355405. This work was supported, in part, by the Cen- 
ter for Functional Nanostructures (CFN) of the Deutsche 
Forschungsgcmcinschaft (DFG) within project C3.9, the 
EXCITiNG Research and Training Network of the Euro- 
pean Union and the NANO QUANTA Network of Excel- 
lence. 



[1] Inhomogeneous electron gas, P. Hohenberg and W. 
Kohn, Phys. Rev. 136, B 864 (1964). 

[2] Self- consistent equations including exchange and corre- 
lation effects, W. Kohn and L.J. Sham, Phys. Rev. 140, 
A 1133 (1965). 

[3] R.M. Dreizler and E.K.U. Gross, Density Functional 



Theory (Springer- Verlag, Berlin, 1990). 
[4] Electron densities in search of Hamiltonians, M. Levy, 

Phys. Rev. A 26, 1200 (1982). 
[5] Density- functional exchange-energy approximation with 

correct asymptotic behavior, A.D. Becke, Phys. Rev. A 

38, 3098 (1988). 



28 



[6] Development of the Colle-Salvetti correlation-energy 
formula into afunctional of the electron density, C. Lee, 
W. Yang, and R.G. Parr, Phys. Rev. B 37, 785 (1988). 

[7] Generalized gradient approximation made simple, [24] 
J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. 
Lett. 77, 3865 (1996); 78, 1396 (1997) (E). 

[8] Density-functional thermochemistry. III. The role of ex- 
act exchange, A.D. Becke, J. Chem. Phys. 98, 5648 
(1993). 

[9] The performance of semi-local and hybrid density func- [25] 
tionals in 3d transition metal chemistry, F. Furche and 
J.P. Perdew, J. Chem. Phys. 124, 044103 (2006). 
[10] see www.turbomole.com 

[11] Time- Dependent Density-Functional Theory, M.A.L. 

Marques and E.K.U. Gross, Annu. Rev. Phys. Chem. [26] 
55, 427 (2004). 

[12] Time- dependent density functional theory in quantum 
chemistry, F. Furche and K. Burke, in Annual Reports 
in Computational Chemistry, Vol 1, ed by D. Spellmeyer [27] 
(Elsevier, Amsterdam, 2005), pg. 19-30. 

[13] Ten topical questions in time- dependent density func- 
tional theory, N. T. Maitra, K. Burke, H. Appel, E.K.U. [28] 
Gross, and R. van Leeuwen, in Reviews in Modern 
Quantum Chemistry: A celebration of the contributions 
of R. G. Parr, ed. K. D. Sen, (World-Scientific, 2001). 

[14] A guided tour of time- dependent density functional the- [29] 
ory, K. Burke and E.K.U. Gross, in Density function- 
als: Theory and applications, ed. D. Joubert (Springer, 
Berlin, 1998). 

[15] Density Functional Theory of Time-Dependent Phe- [30] 
nomena, E.K.U. Gross, J.F. Dobson, and M. Petersilka, 
Topics in Current Chemisty, 181, 81 (1996). 

[16] Density- functional theory for time- dependent systems, 

E. Runge and E.K.U. Gross, Phys. Rev. Lett. 52, 997 [31] 
(1984). 

[17] Dichtefunktionalmethoden fur elektronisch angeregte 
Molekule. Theorie, Implementierung und Anwendung, 

F. Furche, Ph.D. thesis, Universitat Karlsruhe, (2002). 
[18] Auxiliary basis sets for main row atoms and transition 

metals and their use to approximate Coulomb potentials, [32] 
K. Eichkorn, F. Weigend, O. Treutler and R. Ahlrichs, 
Theor. Chem. Acc. 97, 119 (1997). 

[19] Absolute configuration of chiral fullerenes and covalent 
derivatives from their calculated circular dichroism spec- 
tra, H. Goto, N. Harada, J. Crassous and F. Diederich, 
J. Chem. Soc, Perkin Trans. 2, 1719 (1998). [33] 

[20] Pentacyanoiron(II) as an electron donor group for non- 
linear optics: Medium-responsive properties and com- 
parisons with related pentaammineruthenium(II) com- 
plexes, B.J. Coe, J.L. Harries, M. Helliwell, L.A. Jones, 
I. Asselberghs, K. Clays, B.S. Brunschwig, J. A. Har- 
ris, J. Garin and J. Orduna, J. Amer. Chem. Soc. 128, [34] 
12192 (2006). 

[21] Studies of the 5 '-substituted phenylisoquinoline-based 
iridium complexes using density functional theory, C.H. 
Yang, W.L. Su, K.H. Fang, S.P. Wang and I.W. Sun, 
Organometallics 25, 4514 (2006). 

[22] Density-functional theory investigation of the geometric, [35] 
energetic, and optical properties of the cobalt(II)tris(2,2 
'-bipyridine) complex in the high-spin and the John- 
Teller active low-spin states, A. Vargas, M. Zerara, E. 
Krausz, A. Hauser and L.M.L. Daku, J. Chem. Theory [36] 
Comput. 2, 1342 (2006). 

[23] Photon-, electron-, and proton-induced isomerization 



behavior of ferrocenylazobenzenes, A. Sakamoto, A. Hi- 
rooka, K. Namiki, M. Kurihara, M. Murata, M. Sugi- 
moto and H. Nishihara, Inorg. Chem. 44, 7547 (2005). 
Computational and spectroscopic studies of 
Re(I) bipyridyl complexes containing 2,6- 
dimethylphenylisocyanide (CNx) ligand, S.R. Stoyanov, 
J.M. Villegas, A.J. Cruz, L.L. Lockyear, J.H. Reiben- 
spies and D.P. Rillema, J. Chem. Theory Comput. 1, 
95 (2005). 

Theoretical study of Ln(III) complexes with polyaza- 
aromatic ligands: Geometries of [LnL(H20)(n)](3+) 
complexes and successes and failures of TD-DFT, F. 
Gutierrez, C. Rabbe, R. Poteau and J.P. Daudey, J. 
Phys. Chem. A 109, 4325 (2005). 

Importance of vibronic effects on the circular dichro- 
ism spectrum of dimethyloxirane, J. Neugebauer, E.J. 
Baerends, M. Nooijen and J. Autschbach, J. Chem. 
Phys. 122, 234305 (2005). 

Vibronic coupling in square planar complexes of palla- 
dium^!) and platinum(II), E. Lanthier, C. Reber and 
T. Carrington, Chem. Phys. 329, 90 (2006). 
X-ray absorption spectroscopy of titanium oxide by time 
dependent density functional calculations, G. Fronzoni, 
R. Francesco, M. Stener and M. Causa, J. Phys. Chem. 
B 110, 9899 (2006). 

First hyperpolarizability of a sesquifulvalene transition 
metal complex by time-dependent density-functional the- 
ory, W. Hieringer and E.J. Baerends, J. Phys. Chem. A 
110, 1014 (2006). 

Theoretical studies of cyclometalated phenylpyrazol 
Ir(III) complex using density functional theory, N.G. 
Park, G.C. Choi, J.E. Lee and Y.S. Kim, Current Appl. 
Phys. 5, 79 

Synthesis, characterization, and photochemical and 
computational investigations of Ru(II) heterocyclic com- 
plexes containing 2,6-dimethylphenylisocyanide (CNx) 
ligand, J.M. Villegas, S.R. Stoyanov, W. Huang, L.L. 
Lockyear, J.H. Reibenspies and D.P.Rillema, Inorg. 
Chem. 43, 6383 (2004). 

Contrasting linear and quadratic nonlinear optical be- 
havior of dipolar pyridinium chromophores with 4~ 
(dimethylamino)phenyl or ruthenium(II) ammine elec- 
tron donor groups, B.J. Coe, J. A. Harris, B.S. Brun- 
schwig, J. Garin, J. Orduna, S.J. Coles and M.B. Hurst- 
house, J. Am. Chem. Soc. 126, 10418 (2004). 
The charge transfer band solvent- dependence 
of [Ru(bpy)(2)(CNx)Cl](+) CNx=2,6- 

dimethylphenylisocyanide: a polarizable continuum 
model/time- dependent density functional theory study, 
S.R. Stoyanov, J.M. Villegas and D.P. Rillema, Inorg. 
Chem. Comm. 7, 838 (2004). 

Time- dependent density functional theory study of 
the absorption spectrum of [Ru(4,4 '-COOH-2,2 - 
bpy)(2)(NCS)(2)] in water solution: influence of the 
pH, 

F. De Angelis, S. Fantacci and A. Selloni, Chem. Phys. 
Lett. 389, 204 (2004). 

On the origin of the optical activity in the d-d transition 
region of tris-bidentate Co(III) and Rh(III) complexes, 
F.E. Jorge, J. Autschbach and T. Ziegler, Inorg. Chem. 
42, 8902 (2003). 

Simulation of x-ray absorption near-edge spectra and x- 
ray fluorescence spectra of optically excited molecules, 
R.K. Pandey and S. Mukamel, J. Chem. Phys. 124, 



29 



094106 (2006). 

[37] Time-dependent density functional theory calculations 
of ligand K edge and metal L edge X-ray absorption of 
a series of oxomolybdenum complexes, G. Fronzoni, A. 
Stener, A. Reduce and P. Decleva, J. Phys. Chem. A 
108, 8467 (2004). 

[38] TD-DFT investigation on the low-lying excited states 
of spiro-bithiophene, S.Y. Yang, Y.H. Kan, G.C. Yang, 
Z.M. Su and L. Zhao, Chem. Phys. Lett. 429, 180 
(2006). 

[39] Structure absorption spectra correlation in a series 
of 2, 6- dimethyl- 4-arylpyrylium salts, N. Manoj, G. 
Ajayakumar, K.R. Gopidas and C.H. Suresh, J. Phys. 
Chem. A 110, 11338 (2006). 

[40] Optical properties of N-succinimidyl bithiophene and the 
effects of the binding to biomolecules: Comparison be- 
tween coupled-cluster and time- dependent density func- 
tional theory calculations and experiments, E. Fabiano, 
F. Delia Sala, G. Barbarella, S. Lattante, M. Anni, G. 
Sotgiu, C. Hattig, R. Cingolani and G. Gigli, J. Phys. 
Chem. B 110, 18651 (2006). 

[41] Time dependent density functional theory study of elec- 
tronic absorption properties of lead(II) complexes with a 
series of hydroxyflavones, C. Lapouge and J. P. Cornard, 
J. Phys. Chem. A 109, 6752 (2005). 

[42] Optical properties and derealization of excess nega- 
tive charges on oligo(phenylenevinylene)s: A quantum 
chemical study, S. Fratiloiu, F.C. Grozema and L.D.A. 
Siebbeles, J. Phys. Chem. B 109, 5644 (2005). 

[43] A first principles calculations and experimental study 
of the ground- and excited-state properties of ladder 
oligo(p-aniline)s, M. Belletete, G. Durocher, S. Hamel, 
M. Cote, S. Wakim and M. Leclerc, J. Chem. Phys. 122, 
104303 (2005). 

[44] Electronic and optical properties of functionalized car- 
bon chains with the localized Hartree-Fock and conven- 
tional Kohn-Sham methods, M. Weimer, W. Hieringer, 
F. Delia Sala and A. Gorling, Chem. Phys. 309, 77 
(2005). 

[45] Structural, electronic, and optical properties of the di- 
indenoperylene molecule from first-principles density- 
functional theory, L.M. Ramaniah and M. Boero, Phys. 
Rev. A 74, 042505 (2006). 

[46] On the electronic properties of dehydrogenated polycyclic 
aromatic hydrocarbons, D.L. Kokkin and T.W. Schmidt, 
J. Phys. Chem. A 110, 6173 (2006). 

[47] Electronic transitions of thiouracils in the gas phase and 
in solutions: Time-dependent density functional theory 
(TD-DFT) study, M.K. Shukla and J. Leszczynski, J. 
Phys. Chem. A 108, 10367 (2004). 

[48] s-tetrazine in aqueous solution: A density functional 
study of hydrogen bonding and electronic excitations, M. 
Odelius, B. Kirchner and J. Hutter, J. Phys. Chem. A 
108, 2044 (2004). 

[49] Observation and interpretation of annulated por- 
phyrins: Studies on the photophysical properties of 
meso-tetraphenylmetalloporphyrins, J.E. Rogers, K.A. 
Nguyen, D.C. Hufnagle, D.G. McLean, W.J. Su, K.M. 
Gossett, A.R. Burke, SA. Vinogradov, R. Pachter and 
PA. Fleitz, J. Phys. Chem. A 107, 11331 (2003). 

[50] DFT and TDDFT study related to electron transfer in 
nonbonded porphine center dot center dot center dot C- 
60 complexes, T.L.J. Toivonen, T.I. Hukka, O. Cramar- 
iuc, T.T. Rantala and H. Lemmetyinen, J. Phys. Chem. 



A 110, 12213 (2006). 
[51] Optical absorption and electron energy loss spectra of 
carbon and boron nitride nanotubes: a first-principles 
approach, A.G. Marinopoulos, L. Wirtz, A. Marini, V. 
Olevano, A. Rubio and L. Reining, Appl. Phys. A 78, 
1157 (2004). 

[52] Real-time Ab initio simulations of excited carrier dy- 
namics in carbon nanotubes, Y. Miyamoto, A. Rubio 
and D. Tomanek, Phys. Rev. Lett. 97, 126104 (2006). 

[53] Photoelectron spectroscopy of C-84 dianions, O.T. 
Ehrler, J.M. Weber, F. Furche and M.M. Kappes, Phys. 
Rev. Lett. 91, 113006 (2003). 

[54] Photoelectron spectroscopy of fullerene dianions C-76(2- 
), 0-18(2-), and C-84(2-), O.T. Ehrler, F. Furche, J.M. 
Weber and M.M. Kappes, J. Chem. Phys. 122, 094321 
(2005). 

[55] Optical and loss spectra of carbon nanotubes: De- 
polarization effects and intertube interactions, A.G. 
Marinopoulos, L. Reining, A. Rubio, and N. Vast, Phys. 
Rev. Lett. 91, 046402 (2003). 

[56] Structures and electronic absorption spectra of a recently 
synthesised class of photodynamic therapy agents, A.D. 
Quartarolo, N. Russo and E. Sicilia, Chem.-Eur. J. 12, 
6797 (2006). 

[57] Are conical intersections responsible for the ultrafast 
processes of adenine, protonated adenine, and the cor- 
responding nucleosides?, S.B. Nielsen and T.I. Soiling, 
ChemPhysChem 6, 1276 (2005). 

[58] Molecular dynamics in electronically excited states using 
time-dependent density functional theory, I. Tavernelli, 
U.F. Rohrig and U. Rothlisberger, Mol. Phys. 103, 963 
(2005). 

[59] A TDDFT study of the optical response of DNA bases, 
base pairs, and their tautomers in the gas phase, A. Tso- 
lakidis and E. Kaxiras, J. Phys. Chem. A 109, 2373 
(2005). 

[60] Calculating absorption shifts for retinal proteins: Com- 
putational challenges, M. Wanko, M. Hoffmann, P. 
Strodet, A. Koslowski, W. Thiel, F. Neese, T. Frauen- 
heim and M. Elstner, J. Phys. Chem. B 109, 3606 
(2005). 

[61] A time- dependent density functional theory investiga- 
tion of the spectroscopic properties of the beta-subunit in 
C-phycocyanin, Y.L. Ren, J. Wan, X. Xu, Q.Y. Zhang 
and G.F. Yang, J. Phys. Chem. B 110, 18665 (2006). 

[62] Linear and nonlinear optical response of aromatic 
amino acids: A time- dependent density functional in- 
vestigation, J. Guthmuller and D. Simon, J. Phys. 
Chem. A 110, 9967 (2006). 

[63] Time dependent density functional theory modeling of 
chiroptical properties of small amino acids in solution, 
M.D. Kundrat and J. Autschbach, J. Phys. Chem. A 
110, 4115 (2006). 

[64] Theoretical study on photophysical and photosensitive 
properties of aloe emodin, L. Shen, H.F. Ji and H.Y. 
Zhang, Theochem-J. Mol. Struct. 758, 221 (2006). 

[65] Electronically excited states of tryptamine and its micro- 
hydrated complex, M. Schmitt, R. Brause, CM. Marian, 
S. Salzmann and W.L. Meerts, J. Chem. Phys. 125, 
124309 (2006). 

[66] Photoinduced processes in protonated tryptamine, H. 
Kang, C. Jouvet, C. Dedonder-Lardeux, S. Martren- 
chard, C. Charriere, G. Gregoire, C. Desfrancois, J. P. 
Schermann, M. Barat and J. A. Fayeton, J. Chem. Phys. 



30 



122, 084307 (2005). 

[67] A TDDFT study of the excited states of DNA bases and [82 
their assemblies, D. Varsano, R. Di Felice, M.A.L. Mar- 
ques and A. Rubio, J. Phys. Chem. B 110, 7129 (2006). 

[68] TD-DFT calculations of the potential energy curves for 

the trans-cis photo-isomerization of protonated Schiff [83 
base of retinal, H. Tachikawa and T. Iyama, J. Pho- 
tochem. Photobiol. B: Biol. 76, 55 (2004) 

[69] Time- dependent density-functional approach for biologi- 
cal chromophores: The case of the green fluorescent pro- [84 
tein, M. A. L. Marques, X. Lopez, D. Varsano, A. Cas- 
tro, and A. Rubio, Phys. Rev. Lett. 90, 258101 (2003). 

[70] Optical absorption of the blue fluorescent protein: A 
first-principles study, X. Lopez, M.A.L. Marques, A. 
Castro and A. Rubio, J. Am. Chem. Soc. 127, 12329 [85" 

(2005) . 

[71] Single- Reference Calculations of Photochemical Po- 
tential Energy Surfaces: Time- Dependent Density- 
Functional Theory Calculations of the Con- and 
Disrotatory Surfaces for the Woodward-Hoffmann CC [86 
Ring-Opening of Oxirane, F. Cordova, L.J. Doriol, 
A. Ipatov, M.E. Casida and A. Vela, to be published 

(2006) . 

[72] Theoretical predictions for occurrence of charge transfer 

complex within the two synthesized bichromophores con- [87 
sidering the role of their spacers in interactions with the 
pi-orbitals of the redox centers, A.K. De and T. Ganguly, 
Chem. Phys. Lett. 428, 213 (2006). 

[73] Three-dimensional non linear optical chromophores [88 
based on metal-to -lig and charge-transfer from ruthe- 
nium(II) or iron(II) centers, B.J. Coe, J. A. Harris, B.S. 
Brunschwig, I. Asselberghs, K. Clays, J. Garin and J. 
Orduna, J. Amer. Chem. Soc. 127, 13399 (2005). 

[74] NLO properties of metallabenzene-based chromophores: 
A time- dependent density functional study, A. Karton, 
M.A. Iron, M.E. van der Boom and J.M.L. Martin, J. [89" 
Phys. Chem. A 190, 5454 (2005). 

[75] A time dependent density functional theory study of 
alpha-84 phycocyanobilin in C-phycocyanin, J. Wan, X. 
Xu, Y.L. Ren and G.F. Yang, J. Phys. Chem. B 109, 
11088 (2005). 

[76] Photoionization cross section and angular distribution [90 

calculations of carbon tetrafluoride, B. Toffoli, M. 

Stener, G. Fronzoni and P. Decleva, J. Chem. Phys. 

124, 214313 (2006). 
[77] Time-dependent density functional theory determina- [91 

tion of the absorption spectra of naphthoquinones, D. 

Jacquemin, J. Preat, V. Wathelet and E.A. Perpete, 

Chem. Phys. 328, 324 (2006). 
[78] A Theoretical Investigation of the Ground and Excited [92 

States of Coumarin 151 and Coumarin 120, R.J. Cave, 

K. Burke, and E.W. Castner Jr., J. Phys. Chem. A 106, 

9294 (2002). 

[79] A TD-DFT study on triplet excited-state properties of 

curcumin and its implications in elucidating the photo [93 
sensitizing mechanisms of the pigment, L. Shen, H.F. Ji 
and H.Y. Zhang, Chem. Phys. Lett. 409, 300 (2005). 

[80] Time- dependent density functional theory investigation 

of the absorption, fluorescence, and phosphorescence [94 
spectra of solvated coumarins, D. Jacquemin, E.A. Per- 
pete, G. Scalmani, M.J. Frisch, X. Assfeld, I. Ciofini 
and C. Adamo, J. Chem. Phys. 125, 164324 (2006). [95] 

[81] TDDFT investigation of the optical properties of cya- 
nine dyes, B. Champagne, M. Guillaume and F. Zutter- 



man, Chem. Phys. Lett. 425, 105 (2006). 
Excitation spectra of nitro-diphenylaniline: Accurate 
time- dependent density functional theory predictions for 
charge-transfer dyes, D. Jacquemin, M. Bouhy and E.A. 
Perpete, J. Chem. Phys. 124, 204321 (2006). 
Spectral properties of bipyridyl ligands by time- 
dependent density functional theory, F. Labat, P.P. 
Laine, I. Ciofini and C. Adamo, Chem. Phys. Lett. 417, 
445 (2006). 

Photoinduced intramolecular charge transfer reaction in 
(E)-3-(4-Methylamino-phenyl)-acrylic acid methyl es- 
ter: A fluorescence study in combination with TDDFT 
calculation, A. Chakraborty, S. Kar, D.N. Nath and N. 
Guchhait, J. Phys. Chem. A 110, 12089 (2006). 
Theoretical analysis of low-lying charge transfer states 
in [Ru(X) (Me)(CO)(2)(Me-DAB)] (X=Gl, I; DAB=1 
4~diaza-l,3-butadiene) complexes by TDDFT and 
CASSCF/CASPT2 methods, N. Ben Amor, S. Zalis and 

C. Daniel, Int. J. Quant. Chem. 106, 2458 (2006). 
Charge-transfer pi pi* excited state in the 7-azaindole 
dimer. A hybrid configuration interactions singles/time- 
dependent density functional theory description, R. 
Gelabert, M. Moreno and J.M. Lluch, J. Phys. Chem. 
A 110, 1145 (2006). 

Hartree-Fock exchange in time dependent density func- 
tional theory: application to charge transfer excitations 
in solvated molecular systems, L. Bernasconi, M. Sprik 
and R.Hutter, Chem. Phys. Lett. 394, 141 (2004). 
Time- dependent density-functional theory in- 
vestigation of the fluorescence behavior as a 
function of alkyl chain size for the 4~(N,N- 
dimethylamino)benzonitrile-like donor- acceptor sys- 
tems 4-(N,N-diethylamino)benzonitrile and 4~(N,N- 
diisopropylamino)benzonitrile, C.J. Jamorski and M.E. 
Casida, J. Phys. Chem. B 108, 7132 (2004). 
Failure of time- dependent density functional 
theory for long-range charge-transfer excited 
states: The zincbacteriochlorin-bacterlochlorin and 
bacteriochlorophyll-spheroidene complexes, A. Dreuw 
and M. Head-Gordon, J. Am. Chem. Soc. 126, 4007 
(2004). 

Photoinduced intramolecular charge transfer in 4~ 
(dimethyl) aminobenzonitrile - A theoretical perspective, 

D. Rappoport and F. Furche, J. Am. Chem. Soc. 126, 
1277 (2004). 

Relationship between long-range charge-transfer excita- 
tion energy error and integer discontinuity in Kohn- 
Sham theory, D.J. Tozer, J. Chem. Phys. 119, 12697 
(2003). 

Rational classification of a series of aromatic donor- 
acceptor systems within the twisting intramolecu- 
lar charge transfer model, a time- dependent density- 
functional theory investigation, C.J. Jamorski and H.P. 
Luthi, J. Chem. Phys. 119, 12852 (2003). 
Electronic transitions in [Re6S8X6j (4~) (X = CI, Br, 
I): Results from time- dependent density functional the- 
ory and solid-state calculations, L.E. Roy and T. Hugh- 
banks, Inorg. Chem. 45, 8273 (2006). 
Size dependence of the static polarizabilities and absorp- 
tion spectra of Ag-n (n—2-8) clusters, J.C. Idrobo, S. 
Ogut and J. Jellinek, Phys. Rev. B 72, 085445 (2005). 
Finite lifetime effects on the polarizability within time- 
dependent density-functional theory, L. Jensen, J. 
Autschbach, and G.C. Schatz, J. Chem. Phys. 122, 



31 



224115 (2005). 

[96] Polarization effects and charge separation in AgCl-water 

clusters, S.S.M.C. Godinho, P.C. do Couto and B.J.C. 

Cabral, J. Chem. Phys. 122, 044316 (2005). 
[97] Core excitations in MgO: a DFT study with cluster mod- 
els, M. Stener, G. Fronzoni and R. De Francesco, Chem. 

Phys. 309, 49 (2005). 
[98] Absorption spectra of small silver clusters Ag-n (n=4, 

6, 8): A TDDFT study, G.F. Zhao, Y. Lei and Z. Zeng, 

Chem. Phys. 327, 261 (2006). 
[99] Excitonic effects and optical properties of passivated 

CdSe clusters, M.L. del Puerto, M.L. Tiago and J.R. 

Chelikowsky, Phys. Rev. Lett. 97, 096401 (2006). 
[100] Optical properties of sila-adamantane nanoclusters from 

density-functional theory, O. Lehtonen and D. Sund- 

holm, Phys. Rev. B 74, 045433 (2006). 
[101] A theoretical investigation of hyperpolamzability for 

small GanAsm (n+m—4-10) clusters, Y.Z. Lan, W.D. 

Cheng, D.S. Wu, J. Shen, S.P. Huang, H. Zhang, Y.J. 

Gong and F.F. Li, J. Chem. Phys. 124, 094302 (2006). 
[102] Photochemistry of AgCl-water clusters: Comparison 

with Cl-water clusters, S.S.M.C. Godinho, P.C. do 

Couto and B.J.C. Cabral, Chem. Phys. Lett. 419, 340 

(2006). 

[103] Electronic excited states of Si(100) and organic 
molecules adsorbed on Si(100), N.A. Besley and A.J. 
Blundy, J. Phys. Chem. B 110, 1701 (2006). 

[104] Quantum chemical modeling of electrochromism of tung- 
sten oxide films, E. Broclawik, A. Gora, P. Liguzinski, P. 
Petelenz and H.A. Witek, J. Chem. Phys. 124, 054709 
(2006). 

[105] Different approaches for the calculation of electronic ex- 
cited states of nonstoichiometric alkali halide clusters: 
The example of Na s F, G. Durand, M.C. Heitz, F. 
Spiegelman, C. Meier, R. Mitric, V. Bonacic-Koutecky 
and J. Pittner, J. Chem. Phys. 121, 9898 (2004). 

[106] Simple DFT model of clusters embedded in rare gas ma- 
trix: Trapping sites and spectroscopic properties of Na 
embedded in Ar, B. Gervais, E. Giglio, E. Jacquet, A. 
Ipatov, P.G. Reinhard and E. Suraud, J. Chem. Phys. 
121, 8466 (2004). 

[107] Real-time study of the adiabatic energy loss in an atomic 
collision with a metal cluster, R. Baer and N. Siam, J. 
Chem. Phys. 121, 6341 (2004). 

[108] Angular distribution of electrons emitted from Na clus- 
ters, A. Pohl, P.G. Reinhard and E. Suraud, Phys. Rev. 
A 70, 023202 (2004). 

[109] Excited states of tetrahedral single-core S129 nanopar- 
ticles, R. Rao, J. Sutin, R. Clegg, E. Gratton, M.H. 
Nayfeh, S. Habbal, A. Tsolakidis and R.M. Martin, 
Phys. Rev. B 69, 205319 (2004). 

[110] Quantum and thermal fluctuation effects on the pho- 
toabsorption spectra of clusters, F. Delia Sala, R. 
Rousseau, A. Gorling and D. Marx, Phys. Rev. Lett. 
92, 183401 (2004). 

[Ill] Enhanced absorption induced by a metallic nanoshell, 
R. Baer, D. Neuhauser and S. Weiss, Nano Lett. 4, 85 
(2004). 

[112] Photoabsorption spectra of Ti8C12 metallocarbo- 
hedrynes: Theoretical spectroscopy within time- 
dependent density functional theory, J.I. Martinez, A. 
Castro, A. Rubio and J. A. Alonso, J. Chem. Phys. 
125, 074311 (2006). 

[113] Calculation of the optical spectrum of the Ti8C12 and 



V8C12 Met-Cars, J.I. Martinez, A. Castro, A. Rubio, 
J.M. Poblet and J.A. Alonso, Chem. Phys. Lett. 398, 
292 (2004). 

[114] Electronic distribution and solvatochromism investiga- 
tion of a model radical (2,2,6,6-tetramethylpiperidine N- 
oxyl : tempo) through TD-DFT calculations including 
PCM solvation, J. Lalevee, X. Allonas and P. Jacques, 
Theochem-J. Mol. Struct. 767, 143 (2006). 

[115] The merits of the frozen- density embedding scheme to 
model solvatochromic shifts, J. Neugebauer, M.J. Louw- 
erse, E.J. Baerends and T.A. Wesolowski, J. Chem. 
Phys. 122, 094115 (2005). 

[116] Solvatochromism of a novel betaine dye derived from 
purine, A. Masternak, G. Wenska, J. Milecki, B. Skalski 
and S. Franzen, J. Phys. Chem. A 109, 759 (2005). 

[117] On the performance of gradient- corrected approxima- 
tion functionals and polarizable continuum model in the 
study of 1,2,3-triazine in water, C. Zazza and A. Grandi 
and L. Bencivenni and M. Aschi, Theochem-J. Mol. 
Struct. 764, 87 (2006). 

[118] Photophysical properties of trans-3-(4~ 

monomethylamino-phenyl)-acrylonitrile: Evidence 
of twisted intramolecular charge transfer (TICT) pro- 
cess, A. Chakraborty, S. Kar and N. Guchhait, Chem. 
Phys. 324, 733 (2006). 

[119] Water solvent effect on the first hyperpolarizability 
of p-nitrophenol and p-nitrophenylphosphate: A time- 
dependent density functional study, J. Guthmuller and 
D. Simon, J. Chem. Phys. 124, 174502 (2006). 

[120] Geometry and solvent dependence of the electronic spec- 
tra of the amide group and consequences for peptide cir- 
cular dichroism, J. Sebek, Z. Kejik and P. Bour, J. Phys. 
Chem. A 110, 4702 (2006). 

[121] Absorption and emission spectra in gas-phase and solu- 
tion using TD-DFT: Formaldehyde and benzene as case 
studies, D. Jacquemin, E.A. Perpete, G. Scalmani, M.J. 
Frisch, M. Ciofini, C. Adamo, Chem. Phys. Lett. 421, 
272 (2006). 

[122] Geometries and properties of excited states in the gas 
phase and in solution: Theory and application of a time- 
dependent density functional theory polarizable contin- 
uum model, G. Scalmani, M.J. Frisch, B. Mennucci, J. 
Tomasi, R. Cammi and V. Barone, J. Chem. Phys. 124, 
094107 (2006). 

[123] Density functional calculation of the electronic ab- 
sorption spectrum of Cu+ and Ag+ aqua ions, L. 
Bernasconi,J. Blumberger, M. Sprik, and R. Vuilleu- 
mier, J. Chem. Phys. 121, 11885 (2004). 

[124] Magnitude of zero-point vibrational corrections to op- 
tical rotation in rigid organic molecules: A time- 
dependent density functional study, B.C. Mort and J. 
Autschbach, J. Phys. Chem. A 109, 8617 (2005). 

[125] Determination of absolute configurations of chiral 
molecules using ab initio time- dependent density func- 
tional theory calculations of optical rotation: How re- 
liable are absolute configurations obtained for molecules 
with small rotations?, P.J. Stephens, D.M. Mccann, J.R. 
Cheeseman and M.J. Frisch, Chirality 17, S52 (2005). 

[126] A theoretical study of the chiroptical properties of 
molecules with isotopically engendered chirality, M. 
Dierksen and S. Grimme, J. Chem. Phys. 124, 174301 
(2006). 

[127] Time-dependent density functional calculations of opti- 
cal rotatory dispersion including resonance wavelengths 



32 



as a potentially useful tool for determining absolute con- 
figurations of chiral molecules, J. Autschbach, L. Jensen, 
G.C. Schatz, Y.C.E. Tse and M. Krykunov, J. Phys. 
Chem. A 110, 2461 (2006). 
[128] On the origin of optical activity in tris-diamine com- 
plexes of Co(III) and Rh(III): A simple model based on 
time-dependent density function theory, F.E. Jorge, J. 
Autschbach and T. Ziegler, J. Amer. Chem. Soc. 127, 
975 (2005). 

[129] Circular dichroism spectrum of [Co(en)(3)J(3+) in wa- 
ter: A discrete solvent reaction field study, L. Jensen, 
M. Swart, P.T. Van Duijnen and J. Autschbach, Int. J. 
Quant. Chem. 106, 2479 (2006). 

[130] Determination of absolute configuration using concerted 
ab initio DFT calculations of electronic circular dichro- 
ism and optical rotation: Bicyclo[3.3.1]nonane diones, 
P.J. Stephens, D.M. McCann, E. Butkus, S. Stoncius, 
J.R. Cheeseman and M.J. Frisch, J. Org. Chem. 69, 
1948 (2004). 

[131] Calculation of the A term of magnetic circular dichro- 
ism based on time dependent- density functional theory 
I. Formulation and implementation, M. Seth, T. Ziegler, 
A. Banerjee, J. Autschbach, S.J. A. van Gisbergen and 
E.J. Baerends, J. Chem. Phys. 120, 10942 (2004). 

[132] Open shells in reduced-density-matrix-functional theory, 
N.N. Lathiotakis, N. Helbig, and E.K.U. Gross, Phys. 
Rev. A 72, 030501 (R) (2005). 

[133] One-pot synthesis of helical aromatics: Stereoselectiv- 
ity, stability against racemization, and assignment of 
absolute configuration assisted by experimental and the- 
oretical circular dichroism, M. Watanabe, H. Suzuki, 
Y. Tanaka, T. Ishida, T. Oshikawa and A. Tori-i, J. 
Org.Chem. 69 7794 (2004). 

[134] Circular- dichroism studies on artemisinin and epi- 
artemisinin and their beta-cyclodextrin complexes in so- 
lution, G. Marconi, S. Monti, F. Manoli, A.D. Esposti 
and A. Guerrini, Helv. Chim. Acta 87, 2368 (2004). 

[135] Field emission and electronic structures of carbon al- 
lotropes, K. Watanabe, M. Araidai and K. Tada, Thin 
Solid Films 464-65, 354 (2004). 

[136] Effects of the crystal structure in the dynamical electron 
density-response of hep transition metals, G. Gurtubay, 
W. Ku, J.M. Pitarke and A.G. Eguiluz, Comp. Mater. 
Sci. 30, 104 (2004). 

[137] Time-dependent density functional theory of excitation 
energies of closed-shell quantum dots, K. Hirose, Y. 
Meir, and N.S. Wingreen, Physica E 22 486 (2004). 

[138] Long-range contribution to the exchange-correlation ker- 
nel of time- dependent density functional theory, S. 
Botti, F. Sottile, N. Vast, V. Olevano, L. Reining, H.C. 
Weissker, A. Rubio, G. Onida, R. Del Sole and R.W. 
Godby, Phys. Rev. B 69, 155112 (2004). 

[139] Electronic and structural effects on the nonlinear optical 
behavior in push-pull TTF/tricarbonyl chromiun arene 
complexes, B. Insuasty, C. Atienza, C. Seoane, N. Mar- 
tin, J. Garin, J. Orduna, R. Alcala and B. Villacampa, 
J. Org. Chem. 69, 6986 (2004). 

[140] Experimental and ab initio calculational studies on 2,3- 
diketo-benzopiperazine, F.F. Jian and P.S. Zhao, J. Mol. 
Str. 705, 133 (2004). 

[141] Calculations of the third-order nonlinear optical re- 
sponses in push-pull chromophores with a time- 
dependent density functional theory, N. Kobko, A. Ma- 
sunov and S. Tretiak, Chem. Phys. Lett. 392, 444 



(2004) . 

[142] Prediction of two-photon absorption properties for 
organic chromophores using time- dependent density- 
functional theory, A.M. Masunov and S. Tretiak, J. 
Phys. Chem. B 108, 899 (2004). 

[143] Prediction of second-order optical nonlinearity of 
trisorganotin-substituted beta-Keggin polyoxotungstate, 
W. Guan, G.C. Yang, L.K. Yan and Z.M. Su, Inorg. 
Chem. 45, 7864 (2006). 

[144] Theoretical study on the second-order nonlinear optical 
properties of asymmetric spirosilabifluorene derivatives, 
G.C. Yang, Z.M. Su and C.S. Qin, J. Phys. Chem. A 
110, 4817 (2006). 

[145] Origin of the photoinduced optical second harmonic gen- 
eration arising in N-phenyl microcrystalline films, M. 
Makowska-Janusik, S. Tkaczyk and I.V. Kityk, J. Phys. 
Chem. B 110, 6492 (2006). 

[146] Calculation of two-photon absorption spectra of donor- 
pi-acceptor compounds in solution using quadratic re- 
sponse time- dependent density functional theory, P.N. 
Day, K.A. Nguyen and R. Pachter, J. Chem. Phys. 125, 
094103 (2006). 

[147] Calculations of two-photon absorption cross-sections of 
stibene and bis(styryl) benzene derivatives by means of 
TDDFT-SOS method, G.C. Yang, C.S. Qin, Z.M. Su 
and S. Dong, Theochem-J. Mol. Struct. 726, 61 (2005). 

[148] Effects of (multi) branching of dipolar chromophores on 
photophysical properties and two-photon absorption, C. 
Katan, F. Terenziani, O. Mongin, M.H.V. Werts, L. 
Porres, T. Pons, J. Mertz, S. Tretiak and M. Blanchard- 
Desce, J. Phys. Chem. A 109, 3024 (2005). 

[149] TDDFT study of one- and two-photon absorption prop- 
erties: Donor-pi- acceptor chromophores, P.N. Day, K.A. 
Nguyen and R. Pachter, J. Phys. Chem. B 109, 1803 

(2005) . 

[150] Theoretical investigation of two-photon absorption al- 
lowed excited states in symmetrically substituted di- 
acetylenes by ab initio molecular-orbital method, K. 
Ohta and K. Kamada, J. Chem. Phys. 124, 124303 

(2006) . 

[151] Time- dependent density functional theory studies of the 
photo switching of the two-photon absorption spectra in 
stilbene, metacyclophenadiene, and diarylethene chro- 
mophores, A.E. Clark, J. Phys. Chem. A 110, 3790 
(2006). 

[152] Effects of conjugation length, electron donor and accep- 
tor strengths on two-photon absorption cross sections of 
asymmetric zinc-porphyrin derivatives, O. Rubio-Pons, 
Y. Luo and H. Agren, J. Chem. Phys. 124, 094310 
(2006). 

[153] Two-photon excitation of substituted enediynes, J.F. 

Kauffman, J.M. Turner, I.V. Alabugin, B. Breiner, S.V. 

Kovalenko, E.A. Badaeva, A. Masunov and S. Tretiak, 

J. Phys. Chem. A 110, 241 (2006). 
[154] A Primer in Density Functional Theory, ed. C. Fiol- 

hais, F. Nogueira, and M. Marques (Springer- Verlag, 

NY, 2003). 

[155] Universal variational functionals of electron densities, 
first order density matrices, and natural spin-orbitals 
and solution of the v-representability problem M. Levy, 
Proc. Natl. Acad. Sci. (U.S.A.) 76, 6062 (1979). 

[156] Accurate exchange-correlation potentials and total- 
energy components for the helium isoelectronic series, 
C. J. Umrigar and X. Gonze, Phys. Rev. A 50, 3827 



33 



(1994). [174 

[157] Density- functional theory for fractional particle number: 
Derivative discontinuities of the energy, J. P. Perdew, 
R.G. Parr, M. Levy, and J.L. Balduz, Jr., Phys. Rev. 
Lett. 49, 1691 (1982). [175 

[158] Orbital- Dependent Functionals for the Exchange- 
Correlation Energy: A Third Generation of Density 
Functionals, E. Engel, in A Primer in Density Func- [176 
tional Theory, ed. C. Fiolhais, F. Nogueira, and M. 
Marques (Springer- Verlag, NY, 2003), pg 56. [177 

[159] Orbital functionals in density functional theory: the 

optimized effective potential method, T. Grabo, T. [178 
Kreibich, S. Kurth, and E.K.U. Gross, in Strong 
Coulomb correlations in electronic structure: Beyond 
the local density approximation, ed. V.I. Anisimov (Gor- 
don and Breach, Tokyo, 1998). [179 

[160] Orbital- and state- dependent functionals in density- 
functional theory, A. Gorling, J. Chem. Phys. 123, 
062203 (2005). 

[161] Ground state of the electron gas by a stochastic method, [180 

D. M. Ceperley and B. J. Alder, Phys. Rev. Lett. 45, 
566 (1980). 

[162] Self -interaction correction to density-functional approx- [181 
imations for many-electron systems, J. P. Perdew and A. 
Zunger, Phys. Rev. B 23, 5048 (1981). 

[163] Accurate spin- dependent electron liquid correlation en- 
ergies for local spin-density calculations - a critical anal- [182 
ysis, S. H. Vosko, L. Wilk, and M. Nusair, Can. J. Phys. 
58, 1200 (1980). 

[164] Pair distribution function and its coupling- constant av- [183 
erage for the spin-polarized electron gas, J. P. Perdew 
and Y. Wang, Phys. Rev. B 46, 12947 (1992); ibid. 56, 
7018 (1997) (E). [184 

[165] Correlation energy of an electron gas with a slowly vary- 
ing high density, S.-K. Ma and K.A. Brueckner, Phys. 
Rev. 165, 18 (1968). [185 

[166] The adiabatic connection method: A non-empirical hy- 
brid, K. Burke, M. Ernzerhof, and J. P. Perdew, Chem. 
Phys. Lett. 265, 115 (1997). [186 

[167] A new mixing of Hartree-Fock and local density- 
functional theories, A.D. Becke, J. Chem. Phys. 98, 
1372 (1993). [187 

[168] Relevance of the slowly-varying electron gas to atoms, 
molecules, and solids, J. P. Perdew, L. A. Constantin, 

E. Sagvolden, and K. Burke, Phys. Rev. Lett. 97,223002 [188 
(2006). 

[169] Rationale for mixing exact exchange with density func- 
tional approximations, J. P. Perdew, M. Ernzerhof, and [189 
K. Burke, J. Chem. Phys. 105, 9982 (1996). 

[170] Unambiguous exchange- correlation energy density, K. [190 
Burke, F.G. Cruz, and K.C. Lam, J. Chem. Phys. 109, 
8161 (1998). 

[171] Ensemble-Hartree-Fock scheme for excited states. The 
optimized effective potential method, N.I. Gidopoulos, 
P.G. Papaconstantinou and E.K.U. Gross, PHYSICA B [191 
318, 328 (2002). 

[172] Time- dependent optimized effective potential in the lin- 
ear response regime, M. Petersilka, U.J. Gossmann, and 
E.K.U. Gross, in Electronic Density Functional Theory: [192 
Recent Progress and New Directions, eds. J.F. Dobson, 
G. Vignale, and M.P. Das (Plenum, NY, 1998). [193 

[173] Time- dependent optimized effective potential, C.A. Ull- 
rich, U.J. Gossman, and E.K.U. Gross, Phys. Rev. Lett. 
74, 872 (1995). 



Efficient localized Hartree-Fock methods as effective 
exact- exchange Kohn-Sham methods for molecules, F. 
Delia Sala and A. Gorling, J. Chem. Phys. 115, 5718 
(2001). 

The asymptotic region of the Kohn-Sham exchange po- 
tential in molecules, F. Delia Salla and A. Gorling, J. 
Chem. Phys. 116, 5374 (2002). 

Basis set selection for molecular calculations, E.R. 
Davidson and D. Feller, Chem. Rev. 86, 681 (1986). 
Improvements on the direct SCF method, M. Haser and 
R. Ahlrichs, J. Comput. Chem. 10, 104 (1989). 
Neutral Thermochemical Data, H.Y. Afeefy, J.F. Lieb- 
man, and S.E. Stein in NIST Chemistry WebBook 
http://webbook.nist.gov, Eds. P.J. Linstrom and W.G. 
Mallard, June (2005). 

E.K.U. Gross and W. Kohn, in Advances in Quan- 
tum Chemistry, Vol. 21: Density Functional Theory of 
Many-Fermion Systems, edited by S. B. Trickey (Aca- 
demic Press, San Diego, 1990). 

Memory in time-dependent density functional theory, 
N.T. Maitra, K. Burke, and C. Woodward, Phys. Rev. 
Letts, 89, 023002 (2002). 

Demonstration of initial-state dependence in time- 
dependent density functional theory, N.T. Maitra and 
K. Burke, Phys. Rev. A 63, 042501 (2001); 64 039901 
(E). 

A Time- dependent spin-density functional theory for the 
dynamic spin susceptibility, K.L. Liu and S.H. Vosko, 
Can. J. Phys. 67, 1015 (1989). 

I. D'Amico and G. Vignale, Exact exchange-correlation 
potential for a time- dependent two electron system, 
Phys. Rev. B 59, 7876 (1999). 

Correlation in time- dependent density functional theory, 
P. Hessler, N.T. Maitra, and K. Burke, J. Chem. Phys. 
117, 72 (2002). 

Exact Time- Dependent Exchange- Correlation Poten- 
tials for Strong-Field Electron Dynamics, M. Lein, and 
S. Kummel, Phys. Rev. Lett. 94, 143003 (2005). 
Propagators for the time-dependent Kohn-Sham equa- 
tions, A.Castro, M.A.L. Marques and A. Rubio, J. 
Chem. Phys. 121, 3425 (2004). 

Mapping from Densities to Potentials in Time- 
Dependent Density-Functional Theory, R. van Leeuwen, 
Phys. Rev. Lett. 82, 3863, (1999). 

Excitation energies from time- dependent density- 
functional theory, M. Petersilka, U.J. Gossmann, and 
E.K.U. Gross, Phys. Rev. Lett. 76, 1212 (1996). 
Exact conditions in TDDFT, K. Burke, Lect. Notes 
Phys 706, 181 (2006). 

Time- dependent density functional response theory of 
molecular systems: Theory, computational methods, 
and functionals, M.E. Casida, in Recent developments 
and applications in density functional theory, ed. J.M. 
Seminario (Elsevier, Amsterdam, 1996). 
Treatment of electronic excitations within the adiabatic 
approximation of time dependent density functional the- 
ory, R. Bauernschmitt, R. Ahlrichs, Chem. Phys. Lett., 
256, 454 (1996). 

An Undulatory Theory of the Mechanics of Atoms and 
Molecules, E. Schrodinger, Phys. Rev. 28, 1049 (1926). 
Ab Initio Excitation Spectra and Collective Electronic 
Response in Atoms and Clusters, I. Vasiliev, S. Ogiit, 
and J.R. Chelikowsky, Phys. Rev. Lett. 82, 1919 (1999). 



34 



Octopus: a first-principles tool for excited electron-ion 
dynamics, M.A.L. Marques, A. Castro, G.F. Bertsch 
and A. Rubio, Comput. Phys. Commun. 151,60 (2003). 
CPMD, J. Hutter et al., Copyright IBM Zurich Re- [213; 
search Laboratory and MPI fur Festkorperforschung 
1995-2001. 

On the density matrix based approach to time- dependent [2141 
density functional theory, F. Furche, J. Chem. Phys. 
114, 5982 (2001). 

A theoretical study of the electronic spectrum of naph- 
thalene, M. Rubio, M. Merchan, E. Ortf and B.O. Roos, 
Chem. Phys. 179, 395 (1994). 

Hohenberg-Kohn density-functional theory as an im- [215 
plicit Poisson equation for density changes from 
summed fragment densities, A.E. DePristo, Phys. Rev. 
A 54, 3863 (1996). [216; 
Gaussian basis sets of quadruple zeta valence quality for 
atoms H-Kr, F. Weigend, F. Furche and R. Ahlrichs, J. 
Chem. Phys. 119, 12753 (2003). 

Electron affinities of the first-row atoms revisited. Sys- [217] 
tematic basis sets and wave functions, R.A. Kendall, 
T.H. Dunning and R.J. Harrison, J. Chem. Phys. 96, 
6796 (1992). 

Fully optimized contracted Gaussian basis sets for atoms [218 
Li to Kr, A. Schafer, H. Horn and R. Ahlrichs, J. Chem. 
Phys. 97, 2571 (1992). 

Self- consistent molecular orbital methods. XII. Further 
extensions of Gaussian-type basis sets for use in molec- [219 
ular orbital studies of organic molecules, W.J. Hehre, 
R. Ditchfield and J.A. Pople, J. Chem. Phys. 56, 2257 
(1972). 

The influence of polarization functions on molecular or- 
bital hydrogenation energies, P.C. Hariharan and J.A. [220 
Pople, Theoret. Chim. Acta 28, 213 (1973). 
Balanced basis sets of split valence, triple zeta valence 
and quadruple zeta valence quality for H to Rn: Design 
an assessment of accuracy, F. Weigend and R. Ahlrichs, 
Phys. Chem. Chem. Phys. 7, 3297 (2005). 

Gaussian basis sets for use in correlated molecular calcu- [221 
lations. I. The atoms boron through neon and hydrogen, 
T.H. Dunning Jr., J. Chem. Phys. 90, 1007 (1989). 
Fitting the Coulomb potential variationally in linear- 
combination- of- atomic- orbitals density-functional cal- [222 
culations, J.W. Mintmire and B.I. Dunlap, Phys. Rev. 
A 25, 88 (1982). 

Auxiliary basis sets to approximate Coulomb potentials, [223 
K. Eichkorn, O. Treutler, H. Ohm, M. Haser and R. 
Ahlrichs, Chem. Phys. Lett. 240, 283 (1995). 
Analytical time- dependent density functional derivative [224 
methods within the RI-J approximation, an approach to 
excited states of large molecules, D. Rappoport and F. 
Furche, J. Chem. Phys. 122, 064105 (2005). [225] 
Adiabatic time- dependent density functional methods for 
excited state properties, F. Furche and R. Ahlrichs, J. 
Chem. Phys. 117, 7433 (2002). 

The second-order approximate coupled cluster singles [226 

and doubles model CC2, O. Christiansen, H. Koch and 

P. J0rgensen, Chem. Phys. Lett. 243, 409 (1995). 

CC2 excitation energy calculations on large molecules 

using the resolution of the identity approximation, C. 

Hattig and F. Weigend, J. Chem. Phys. 113, 5154 [227 

(2000). 

Transition moments and excited state first-order proper- 

[228; 



ties in the second-order coupled cluster model CC2 using 
the resolution of the identity approximation, C. Hattig 
and A. Kohn, J. Chem. Phys. 117, 6939 (2002). 
A TDDFT study of the lowest excitation energies of 
polycyclic aromatic hydrocarbons, M. Parac and S. 
Grimme, Chem. Phys. 292, 11 (2003). 
On the required shape corrections to the local den- 
sity and generalized gradient approximations to the 
Kohn-Sham potentials for molecular response calcula- 
tions of (hyper)polarizabilities and excitation energies, 
M. Griming, O.V. Gritsenko, S.J. A. van Gisbergen, and 
E.J. Baerends, J. Chem. Phys. 116, 9591 (2002). 
Density- functional theory calculations with correct long- 
range potentials,Q. Wu, P.W. Ayers, and W. Yang, J. 
Chem. Phys. 119, 2978 (2003). 

Excitation energies of molecules by time- dependent den- 
sity functional theory based on effective exact exchange 
Kohn-Sham-potentials, F. Delia Sala and A. Gorling, 
Int. J. Quantum Chem. 91, 131 (2003). 
Medium-size polarized basis-sets for high-level-correlated 
calculations of molecular electric properties .2. 2nd-row 
atoms - Si through CI, A.J. Sadlej, Theor. Chim. Acta 
79, 123 (1991). 

Molecular excitation energies from time- dependent den- 
sity functional theory, T. Grabo, M. Petersilka, and 
E.K.U. Gross, in J. Mol. Structure (Theochem), 501, 
353 (2000). 

Dynamic polarizabilities and excitation spectra from 
a molecular implementation of time- dependent density 
functional response theory: N2 as a case study, C. 
Jamorski, M.E. Casida, D.R. Salahub, J. Chem. Phys. 
104, 5134 (1996). 

Molecular excitation energies to high-lying bound states 
from time- dependent density-functional response theory: 
Characterization and correction of the time- dependent 
local density approximation ionization threshold M.E. 
Casida, C. Jamorski, K.C. Casida, and D.R. Salahub, 
J. Chem. Phys. 108, 4439 (1998). 

Comparison between equation of motion and polar- 
ization propagator calculations, J. Oddershede, N.E. 
Griiner and G.H.F. Diercksen, Chem. Phys. 97, 303 
(1985). 

K.P. Huber and G. Herzberg, Molecular Spectra 
and Molecular Structure IV: Constants of Diatomic 
Molecules, Van Nostrand Reinhold, New York, 1979. 
Relationship of Kohn-Sham eigenvalues to excitation 
energies, A. Savin, C.J. Umrigar, X. Gonze, Chem. 
Phys. Lett. 288, 391 (1998). 

The development of new exchange-correlation function- 
al, D.J. Tozer and N.C. Handy, J. Chem. Phys. 108, 
2545 (1998). 

A direct optimization method for calculating density 
functionals and exchange- correlation potentials from 
electron densities, Q. Wu and W. Yang, J. Chem. Phys. 
118, 2498 (2003). 

Asymptotic correction approach to improving approxi- 
mate exchange-correlation potentials: Time- dependent 
density-functional theory calculations of molecular ex- 
citation spectra, M.E. Casida, D.R. Salahub, J. Chem. 
Phys. 113, 8918 (2000). 

Exchange-correlation potential with correct asymptotic 
behavior, R. van Leeuwen and E.J. Baerends, Phys. Rev. 
A. 49, 2421 (1994). 

J. B. Krieger, Y. Li, and G. J. Iafrate, in Density Func- 



35 



tional Theory, eds. R. Dreizler and E. K. U. Gross, [248 
NATO ASI Series B (Plenum, New York, 1995). 
Excitation energies from time- dependent density func- 
tional theory using exact and approximate functionals, [249 
M. Petersilka, E.K.U. Gross, and K. Burke, Int. J. 
Quantum Chem. 80, 534 (2000). [250 
First-principles density-functional calculations for opti- 
cal spectra of clusters and nanocrystals, I. Vasiliev, S. 
Ogiit, and J.R. Chelikowsky, Phys. Rev. B 65, 115416 [251 
(2002). 

Excitations in Time- Dependent Density-Functional 
Theory, H. Appel, E.K.U. Gross, and K. Burke, Phys. [252 
Rev. Lett. 90, 043005 (2003). 

Double-Pole Approximation in Time- Dependent Density [253 

Functional Theory H. Appel, E.K.U. Gross, K. Burke, 

Int. J. Quant. Chem. 106, 2840-2847 (2006). 

Measuring the kernel of time- dependent density func- [254 

tional theory with X-ray absorption spectroscopy of 3d 

transition metals, A. Scherz, E.K.U. Gross, H. Appel, 

C. Sorg, K. Baberschke, H. Wende, and K. Burke, Phys. [255 

Rev. Lett. 95, 253006 (2005) 

Comparison of approximate and exact density function- 
als: A quantum monte carlo study, C. J. Umrigar and 
X. Gonze, in High Performance Computing and its Ap- [256 
plication to the Physical Sciences, Proceedings of the 
Mardi Gras 1993 Conference, edited by D. A. Browne [257 
et al. (World Scientific, Singapore, 1993). 
The quantum defect method, F. Ham, Solid State Phys. 
1, 127 (1955). 

The quantum defect method, M. J. Seaton, Mon. Not. [258 
R. Astron. Soc. 118, 504 (1958). 

Optimized effective atomic central potential, J.D. Tal- 
man and W.F. Shadwick, Phys. Rev. A 14, 36 (1976). [259 
New Kohn-Sham method for molecules based on an ex- 
change charge..., A. Gorling, Phys. Rev. Letts. 83, 5459 
(1999). [260 
Exact exchange treatment for molecules in finite-basis- 
set Kohn-Sham theory, S. Ivanov, S. Hirata, and R.J. 
Bartlett, Phys. Rev. Letts. 83, 5455 (1999). [261 
Shape corrections to exchange- correlation potentials by 
gradient-regulated seamless connection of model poten- 
tials for inner and outer region, M. Gruning, O.V. Grit- [262 
senko, S.J. A. van Gisbergen, and E.J. Baerends, J. 
Chem. Phys. 114, 652 (2001). 

The quantum defect: the true measure of TDDFT re- [263 
suits for atoms, M. van Faassen and K. Burke, J. Chem. 
Phys. 124, 094102 (2006). 

The quantum defect as a powerful tool for studying Ry- [264 

dberg transition energies from density functional theory, 

M. van Faassen, Int. J. Quant. Chem. 106, 3235 (2006). 

A new challenge for time-dependent density functional [265 

theory, M. van Faassen and K. Burke, Chem. Phys. Lett. 

431, 410 (2006). 

The role of exchange and correlation in time-dependent [266 
density- functional theory for photoionization, M. 
Stener, P. Decleva, and A. Gorling, J. Chem. Phys., 
114, 7816 (2001). 

Accurate oscillator strengths for neutral helium, A. [267 
Kono and S. Hattori, Phys. Rev. A 29, 2981 (1984). 
Precision measurements of the absolute photoionization [268 
cross sections of He, J.A.R. Samson, Z.X. He, L. Yin 
and G.N. Haddad, J. Phys. B 27, 887 (1994). 
Theoretical Atomic Physics, H. Friedrich, 2nd ed. [269 
(Springer- Ver lag, 1998). 



Rydberg transition frequencies from the Local Density 
Approximation, A. Wa sserman and K. Burke, Phys. 
Rev. Lett. 95, 163006 (2005). 

Handbook of Mathematical Functions, eds. M. 
Abramowitz and LA. Stegun (Dover, New York, 1972). 
Accurate Rydberg transitions from LDA potentials, A. 
Wasserman, N.T. Maitra, and K. Burke, Phys. Rev. 
Lett. 91, 263001 (2003). 

Electron-molecule scattering from time-dependent den- 
sity functional theory, A. Wasserman, N.T. Maitra, and 
K. Burke, J. Chem. Phys. 122, 133103(2005). 
Scattering amplitudes from TDDFT, A. Wasserman and 
K. Burke, Lect. Notes Phys 706, 493 (2006). 
Bound-free correlation in electron scattering by atoms 
and molecules, R.K. Nesbet, Phys. Rev. A 62, 
040701(R) (2000). 

Electron-molecule scattering calculations in a 3D fi- 
nite element R-matrix approach, S. Tonzani and C. H. 
Greene, J. Chem. Phys. 122, 014111 (2005). 
Low-energy electron scattering from DNA and RNA 
bases: Shape resonances and radiation damage, S. Ton- 
zani and C. H. Greene, J. Chem. Phys. 124, 054312 
(2006). 

On the absorption spectrum of noble gases at the arc 
spectrum limit, U. Fano, cond-mat/0502210, (2005). 
Excitation energies from time- dependent density- 
functional theory beyond the adiabatic approximation, 
C.A. Ullrich and K. Burke, J. Chem. Phys. 121, 28 

(2004) . 

Nonadiabatic electron dynamics in time- dependent 
density-functional theory, C.A. Ullrich and I.V. Tokatly, 
Phys. Rev. B 73, 235102 (2006). 

Time-dependent density-functional theory beyond the 
adiabatic approximation: insights from a two-electron 
model system, C. A. Ullrich, cond-mat/0610341 (2006). 
Time- dependent density functional theory beyond the 
adiabatic approximation, G. Vignale, Int. J. Mod. Phys. 
B 15, 1714 (2001). 

Time- dependent density-functional theory beyond the 
local-density approximation, J.M. Tao and G. Vignale, 
Phys. Rev. Lett. 97, 036403 (2006). 
Time-dependent exchange-correlation current density 
functionals with memory, Y. Kurzweil and R. Baer, J. 
Chem. Phys. 121, 8731 (2004). 

Generic Galilean-invariant exchange-correlation func- 
tionals with quantum memory, Y. Kurzweil and R. Baer, 
Phys. Rev. B 72, 035106 (2005). 

Quantum memory effects in the dynamics of electrons 
in gold clusters, Y. Kurzweil and R. Baer, Phys. Rev. B 
73, 075413 (2006). 

Quantum many-body dynamics in a Lagrangian frame: 

I. Equations of motion and conservation laws, I.V. 
Tokatly, Phys. Rev. B 71, 165104 (2005). 

Quantum many-body dynamics in a Lagrangian frame: 

II. Geometric formulation of time- dependent density 
functional theory, I.V. Tokatly, Phys. Rev. B 71, 165105 

(2005) . 

Current Density Functional Theory G. Vignale, Lect. 
Notes Phys 706, 75 (2006). 

Current-Dependent Exchange- Correlation Potential for 
Dynamical Linear Response Theory, G. Vignale and W. 
Kohn, Phys. Rev. Lett. 77, 2037 (1996). 
Time-dependent density functional theory beyond the 
adiabatic local density approximation, G. Vignale, C.A. 



36 



Ullrich, and S. Conti, Phys. Rev. Lett. 79, 4878 (1997). [287 

[270] Double excitations in time- dependent density functional 
theory linear response, N.T. Maitra, F. Zhang, R.J. 
Cave and K. Burke, J. Chem. Phys. 120, 5932 (2004). 

[271] Approximate Functionals from Many-Body Perturbation 
Theory, A. Marini, R. Del Sole and A. Rubio, Lect. 
Notes Phys 706, 161 (2006). 

[272] Does Density Functional Theory Contribute to the Un- [289 
derstanding of Excited States of Unsaturated Organic 
Compounds?, D.J. Tozer, R.D. Amos, N.C. Handy, 
B.O. Roos, and L. Serrano-Andres, Mol. Phys. 97, 859 
(1999). [290 

[273] On the determination of excitation energies using den- 
sity functional theory, D.J. Tozer and N.C. Handy, Phys. 
Chem. Chem. Phys. 2, 2117 (2000). [291 

[274] Adiabatic connection for near degenerate excited states, 
F. Zhang and K. Burke, Phys. Rev. A 69, 052510 (2004). 

[275] Propagator corrections to adiabatic time-dependent [292 
density-functional theory linear response theory, M.E. 
Casida, J. Chem. Phys. 122, 054111 (2005). 

[276] A dressed TDDFT treatment of the 2 1 A g states of buta- [293 
diene and hexatriene, R.J. Cave, F. Zhang, N.T. Maitra, 
and K. Burke, Chem. Phys. Lett. 389, 39 (2004) 

[277] Resonant nonlinear polarizabilities in the time- [294 
dependent density functional (TDDFT) theory, S. 
Tretiak and V. Chernyak, J. Chem. Phys. 119, 8809 
(2003). [295 

[278] Ultranonlocality in time- dependent current- density- 
functional theory: application to conjugated polymers, 
M. van Faassen, P.L. de Boeij, R. van Leeuwen, J. A. 
Berger, and J.G. Snijders, Phys. Rev. Lett. 88, 186401 
(2002). 

[279] Application of time- dependent current- density- [296 
functional theory to nonlocal exchange-correlation 
effects in polymers, M. van Faassen, P.L. de Boeij, R. 
van Leeuwen, J. A. Berger, and J.G. Snijders, J. Chem. [297 
Phys. 118, 1044 (2003). 

[280] Time- dependent current-density-functional theory ap- 
plied to atoms and molecules, M. van Faassen, Int. J. [298 
Mod. Phys. B 20, 3419 (2006). 

[281] Electric field dependence of the exchange-correlation 
potential in molecular chains, S.J. A. van Gisbergen, 
P.R.T. Schipper, O.V. Gritsenko, E.J. Baerends, J.G. [299 
Snijders, B. Champagne, B. Kirtman, Phys. Rev. Lett., 
83, 694 (1999). 

[282] Electrical Response of Molecular Chains from Density [300 
Functional Theory, S. Kummel, L. Kronik, and J. P. 
Perdew, Phys. Rev. Lett. 93, 213002 (2004). 

[283] Optical excitations of Si by time- dependent density func- 
tional theory based on exact- exchange Kohn-Sham band [301 
structure, Y. H. Kim, M. Stadele, and A. Gorling, Int. 
J. Quantum Chem. 91, 257 (2003). 

[284] Efficient real-space approach to time- dependent density [302 
functional theory for the dielectric response of nonmetal- 
lic crystals, F. Kootstra, P.L. Boeij, and J.G. Snijders, 
J. Chem. Phys. 112, 6517 (2000). [303 

[285] Application of time- dependent density-functional theory 
to the dielectric function of various nonmetallic crystals, 
F. Kootstra, P. L. De boeij, and J. G. Snijders, Phys. 
Rev. B, 62, 7071 (2000). [304 

[286] Current density functional theory for optical spectra: A 
polarization functional, P.L. de Boeij, F. Kootstra, J. A. 
Berger, R. van Leeuwen, and J.G. Snijders, J. Chem. [305 
Phys. 115, 1995 (2001). 



Electronic excitations: density-functional versus many- 
body Green's- function approaches, G. Onida, L. Reining, 
and A. Rubio, Rev. Mod. Phys. 74, 601 (2002). 
Parameter- free calculation of response functions in 
time- dependent density-functional theory, F. Sottile, V. 
Olevano, and L. Reining, Phys. Rev. Lett. 91, 056402 
(2003). 

Bound Excitons in Time- Dependent Density- Functional 
Theory: Optical and Energy-Loss Spectra, Andrea 
Marini, Rodolfo Del Sole, and Angel Rubio, Phys. Rev. 
Lett. 91, 256402 (2003) 

Electron linewidths of wide-gap insulators: Excitonic ef- 
fects in LiF, A. Marini and A. Rubio, Phys. Rev. B 70, 
081103 (2004). 

Time-dependent density functional theory, edited by 
M. Marques, C.A. Ullrich, F. Noguiera, A. Rubio, K. 
Burke, and E.K.U. Gross (Springer, Heidelberg, 2006). 
Undoing static correlation: long-range charge trans- 
fer in time-dependent density functional theory, N. T. 
Maitra, J. Chem. Phys. 122, 234104 (2005). 
Long-range excitations in time- dependent density func- 
tional theory, N.T. Maitra, and D.G. Tempel, cond- 
mat/0510588 (2005). 

Failure of time- dependent density functional methods for 
excitations in spatially separated systems, W. Hieringer, 
and A. Gorling, Chem. Phys. Lett. 419, 557 (2006). 
Charge-transfer correction for improved time- dependent 
local density approximation excited-state potential en- 
ergy curves: Analysis within the two-level model with il- 
lustration for H-2 and LiH, M.E. Casida, F. Gutierrez, 
J.G. Guan, F.X. Gadea, D. Salahub, and J. P. Daudey, 
J. Chem. Phys. 113, 7062 (2000). 

Fluctuation-dissipation theorem density-functional the- 
ory, F. Furche and T. Van Voorhis, J. Chem. Phys. 122, 
164106 (2005). 

Density functional approximation for the correlation en- 
ergy of the inhomogeneous gas, J. P. Perdew, Phys. Rev. 
B 33, 8822 (1986); 34, 7406 (1986) (E). 
Electron correlation energies from scaled exchange- 
correlation kernels: Importance of spatial vs. temporal 
nonlocality, M. Lein, E.K.U. Gross, and J. P. Perdew, 
Phys. Rev. B 61, 13431 (2000). 

van der Waals Interactions in Density-Functional The- 
ory, Y. Andersson, D.C. Langreth, and B.I. Lunqvist, 
Phys. Rev. Lett. 76, 102 (1996). 

Describing static correlation in bond dissociation by 
Kohn-Sham density fun ctional theory, M. Fuchs, Y.- 
M. Niquet, X. Gonze, and K. Burke, J. Chem. Phys. 
122, 094116 

van der Waals density functional for general geometries, 
M. Dion, H. Rydberg, E. Schrder, D.C. Langreth, and 
B.I. Lundqvist, Phys. Rev. Lett. 92, 246401 (2004). 
The exchange-correlation energy of a metallic surface, 
D.C. Langreth and J. P. Perdew, Solid State Commun. 
17, 1425 (1975). 

Exchange and correlation in atoms, molecules, and 
solids by the spin- density-functional formalism, O. Gun- 
narsson and B.I. Lundqvist, Phys. Rev. B 13, 4274 
(1976). 

Density Functional for van der Waals Forces at Sur- 
faces, Y. Andersson, B.I. Lunqvist, and D.C. Langreth, 
Phys. Rev. Lett. 77, 2029 (1996). 

Density functional results for isotropic and anisotropic 



37 



multipole polarizabilities and C%, CV, and Cg Van der 
Waals dispersion coefficients for molecules, V.P. Osinga, 
S.J. A. van Gisbergen, J.G. Snijders, and E.J. Baerends, [324] 
J. Chem. Phys. 106, 5091 (1997). [325] 

[306] Van der Waals energies in density functional theory, W. 
Kohn, Y. Meir, and D.E. Makarov, Phys. Rev. Lett., 80, 
4153 (1998). [326] 

[307] Van der Waals interactions in density functional the- 
ory, Y. Andersson, E. Hult, H. Rydberg, P. Apell, 
B.I. Lundqvist, and D.C Langreth, in Electronic Den- [327] 
sity Functional Theory: Recent Progress and New Di- 
rections, eds. J.F. Dobson, G. Vignale, and M.P. Das 
(Plenum, NY, 1997). 

[308] Tractable nonlocal correlation density functionals for 

fiat surfaces and slabs, H. Ryberg, B.I. Lundqvist, D.C. [328] 
Langreth, and M. Dion, Phys. Rev. B 62, 6997 (2000). 

[309] Van der Waals Density Functional for Layered Struc- 
tures, H. Rydberg, M. Dion, N. Jacobson, E. Schroder, 
P. Hyldgaard, S.I. Simak, D.C. Langreth, and B.I. 
Lundqvist, Phys. Rev. Lett, (to appear). [329] 

[310] Application of van der Waals Density Functional to 
an Extended System: Adsorption of Benzene and 
Naphthalene on Graphite, S.D. Chakarova Kack, E. 
Schroder, B.I. Lundqvist, and D.C. Langreth, cond- 
mat/0601161vl (2006). [330] 

[311] Binding energies in benzene dimers: Nonlocal density 
functional calculations, A. Puzder, M. Dion, and D.C. 
Langreth, J. Chem. Phys. 124, 164105 

[312] Interaction energies of mono substituted benzene dimers 

via nonlocal density functional theory, T. Thonhauser, [331] 
A. Puzder, and D.C. Langreth, J. Chem. Phys. 124, 
164106 

[313] Dispersion Energy from Density-Functional Theory De- 
scription of Monomers, A.J. Misquitta, B. Jeziorski, [332] 
and K. Szalewicz, Phys. Rev. Lett. 91, 033201 (2003). 

[314] Molecular tests of the random phase approximation to 
the exchange-correlation energy functional F. Furche, 
Phys. Rev. B 64, 195120 (2001). [333] 

[315] Accurate density functionals: Approaches using the [334] 
adiabatic- connection fluctuation- dissipation theorem, 
M. Fuchs and X. Gonze, Phys. Rev. B 65, 235109 
(2002). 

[316] Double ionization of helium at 390 nm, J. Parker, [335] 
L.R. Moore, D. Dundas and K.T. Taylor, J. Phys. B: 
At. Mol. Opt. Phys. 33, L691-L698 (2000). 

[317] Time-Dependent Variational Approach to Molecules in 

Strong Laser Fields, T. Kreibich, R. van Leeuwen, and [336] 
E.K.U. Gross, Chem. Phys. 304, 183 (2004). 

[318] Towards time- dependent density-functional theory for 
molecules in strong laser pulses, T. Kreibich, N.I. Gi- 
dopoulos, R. van Leeuwen, E.K.U. Gross, Progress in [337] 
Theoretical Chemistry and Physics 14, 69 (2003). 

[319] Strong-Field double Ionization of Helium: A Density- 
Functional perspective, M. Petersilka and E.K.U. Gross, [338] 
Laser Physics 9,1 (1999). 

[320] The Effect of the Electron- Electron Interaction in 
Above-Threshold Double Ionization, M. Lein, E.K.U. 
Gross, and V. Engel, Laser Physics 12, 487 (2002). [339] 

[321] Strong-field ionization dynamics of a model H2 molecule, 

M. Lein, T. Kreibich, E.K.U Gross, and V. Engel, Phys. [340] 
Rev. A 65, 033403 (2002). 

[322] Phase-space analysis of double ionization, M. Lien, V. 

Engel and E.K.U. Gross, Optics Express 8, 411 (2001). [341] 

[323] Intense-Field Double Ionization of Helium: Identifying 



the Mechanism, M. Lein, E.K.U. Gross, and V. Engel, 

Phys. Rev. Lett. 85, 4707 (2000). 

http:/ / www. physik.fu- berlin.de/ ~ag-gross/tdelf 

Time- dependent electron localization function, T. Bur- 

nus, M.A.L. Marques, E.K.U. Gross, Phys. Rev. A 71, 

010501 (R) (2005). 

Time- dependent electron localization functions for cou- 
pled nuclear- electronic motion, M. Erdmann, E.K.U. 
Gross, V. Engel, J. Chem. Phys. 121, 9666 (2004). 
Time- dependent density-functional approach to atoms 
in strong laser pulses, C.A. Ullrich, S. Erhard and 
E.K.U. Gross, in: Super Intense Laser Atom Physics IV, 
H.G. Muller and M.V. Fedorov, ed(s), (Kluwer, 1996), 
p 267 - 284. 

High harmonic generation in hydrogen and helium 
atoms subject to one- and two-color laser pulses, S. 
Erhard and E.K.U. Gross, in: Multiphoton Processes 
1996, P. Lambropoulus and H. Walther, ed(s), (IOP, 
1997), p 37. 

Excited states dynamics in time- dependent density func- 
tional theory - High-field molecular dissociation and har- 
monic generation, A. Castro, M.A.L. Marques, J. A. 
Alonso, G.F. Bertsch and A. Rubio, Eur. Phys. J. D 
28, 211 (2004). 

Self -interaction- free time- dependent density-functional 
theory for molecular processes in strong fields: high- 
order harmonic generation of Hi in intense laser fields 
- art. no. 023+11, X. Chu and S. I. Chu, Phys. Rev. A, 
6302, 3411 (2001). 

Ionization and High-order harmonic generation in 
aligned benzene by a short intense circularly polarized 
laser pulse, R. Baer, D. Neuhauser, P. R.Zdanska and 
N. Moiseyev, Phys. Rev. A 68,043406 (2003). 
Monitoring molecular dynamics using coherent electrons 
from high harmonic generation, N.L. Wagner, A. Wiiest, 
LP. Christov, T. Popmintchev, X. Zhou, M.M. Murnane 
and H.C. Kapteyn, PNAS 103, 13279 (2006). 
see http://online.itp.ucsb.edu/online/atto06/moletomo/ 
Precision measurement of stong field double ionization 
of helium, B. Walker, B. Sheehy, L.F. DiMauro, P. Agos- 
tini, K.J. Schafer, and K.C. Kulander, Phys. Rev. Lett. 
73, 1227 (1994). 

Observation of nonsequential double ionization of He- 
lium with optical tunnelling, D.N. Fittinghoff, P.R. 
Bolton, B. Chang, and K.C. Kulander, Phys. Rev. Letts. 
69, 2642 (1992). 

Wavelength dependence of nonsequential double ioniza- 
tion in He , K. Kondo, A. Sagisaka, T. Tamida, Y. 
Nabekawa and S. Watanabe, Phys. Rev. A 48, R2531 
(1993). 

Towards Single-Particle spectroscopy of Small Metal 
Clusters, A. Pohl, P.-G. Reinhard and E. Suraud, Phys. 
Rev. Lett 84,5090 (2000). 

Asymmetry of above-threshold ionization of metal clus- 
ters in two-color laser fields: A time- dependent density- 
functional study, H.S. Nguyen, A.D. Bandrauk and C.A. 
Ullrich, Phys. Rev. A 69, 063415 (2004). 
Optical Control of Molecular Dynamics, S.A. Rice and 
M. Zhao (John Wiley and Sons, 2000). 
Tailoring laser pulses with spectral and fiuence con- 
straints using optimal control theory, J. Werschnik, and 
E.K.U. Gross, J. Opt. B: Quantum Semiclass. Opt. 7, 
Laser-induced control of (multichannel) intracluster re- 
actions - The slowest is always the easiest to take, A.G. 



38 



Urena, K. Gasmi, S. Skowronek, A. Rubio and P.M. 
Echenique, Eur. Phys. J. D 28, 193 (2004). [353] 
Coherent control of intersubband optical bistability in 
quantum wells, FLO. Wijewardane and C.A. Ullrich, 
Appl. Phys. Lett. 84, 3984 (2004). [354] 
Theory of the linewidth of intersubband plasmons in 
quantum wells, C.A. Ullrich and G. Vignale, Phys. Rev. 
Lett. 87, 037402 (2001). [355] 
Collective charge-density excitations of noncircular 
quantum dots in a magnetic field, C. A. Ullrich and G. 
Vignale, Phys. Rev. B, 61, 2729 (2000). 

Theoretical studies of molecular scale near-field electron [356] 
dynamics, R. Baer and D. Neuhauser, J. Chem. Phys. 
125, 074709 (2006). 

Time-dependent density-functional studies of the D-2 
Coulomb explosion, E. Livshits and R. Baer, J. Phys. [357] 
Chem. A 110, 8443 (2006). 

Rotational aspects of short-pulse population transfer in 
diatomic molecules, Y.B. Band, S. Kalhish and R. Baer, [358] 
Chem. Phys. Lett. 392, 23 (2004). 

Electron transport in molecular wire junctions, A. 
Nitzan and M.A. Ratner, Science 300, 1384 (2003). 
Density functional method for nonequilibrium electron [359] 
transport M. Brandbyge, J.-L. Mozos, P. Ordejon, J. 
Taylor, and K. Stokbro, Phys. Rev. B 65, 165401 (2002). 
Time- Dependent Partition- Free Approach in Resonant 
Tunneling Systems, G. Stefanucci and C- O. Almbladh, [360] 
Phys. Rev. B 69, 195318 (2004). 

Transport in nanoscale systems: the microcanonical 

versus grand-canonical picture, M. Di Ventra, T.N. [361] 

Todorov, J Phys Cond Matt 16, 8025 (2004). 

Kinetic theory of quantum transport at the nanoscale, R. 



Gebauer and R. Car, Phys. Rev. B 70, 125324 (2004). 
Density Functional Theory of the Electrical Conductiv- 
ity of Molecular Devices, K. Burke, Roberto Car, and 
Ralph Gebauer, Phys. Rev. Lett. 94, 146803 (2005). 
The conductance of molecular wires and DFT based 
transport calculations, F. Evers, F. Weigend, M. Koen- 
topp, Phys Rev. B 69, 25411 (2004). 
Dynamical Corrections to the DFT-LDA Electron Con- 
ductance in Nanoscale Systems, N. Sai, M. Zwolak, G. 
Vignale, and M. Di Ventra, Phys. Rev. Lett. 94, 186810 

(2005) . 

The Dramatic Role of the Exchange-Correlation Poten- 
tial in ab initio Electron Transport Calculations, S-H 
Ke, H.U. Baranger and W. Yang, cond-mat/0609637 

(2006) . 

Self -interaction errors in density functional calculations 
of electromctransport, C. Toher, A. Filippetti, S. San- 
vito, and K. Burke, Phys. Rev. Lett. 95, 146402 (2005) 
Zero-bias molecular electronics: Exchange-correlation 
corrections to Landauer's formula, M. Koentopp, K. 
Burke, and F. Evers, Phys. Rev. B Rapid Comm., 73, 
121403 (2006). 

Time- dependent quantum transport: A practical scheme 
using density functional theory, S. Kurth, G. Stefanucci, 
CO. Almbladh, A. Rubio and E.K.U. Gross, Phys. Rev. 
B 72, 035308 (2005). 

Kohn-Sham master equation approach to transport 
through single molecules, R. Gebauer, K. Burke, and 
R. Car, Lect. Notes Phys 706, 463 (2006). 
Electron transport with dissipation: a quantum kinetic 
approach R. Gebauer and R. Car, Int. J. Quantum 
Chem. 101, SI 564, (2005). 



